\: 


Wee 


Gx aapais 
UNIERSITATIS 
AARERTAEDSTS 


Digitized by the Internet Archive 
In 2022 with funding from 
University of Alberta Library 


https://archive.org/details/Yung 1984 


ap i 


at an iq 


at Se ¢ lo a ly 7 » 
\ 7 . oat) : Tol a 
( ‘16 


bis ae ae 
ars 


ae le am 


tot a i 
cot 


Yirery wey ele ae |? 
i 2h SOLARA TROR | aa 1 ; 
Se yh wAR sever itu ahabieiel il netnCE 
Tie. FALG ra Rite 


p 
w, 


fie 


Me eraby Granted £6 Teg we ie 
eroovucw #ieG.« “Be oe WV. Ae 
Par € cueh HOples Raye. VF aot ¥) 


ws en . Y 1OQORSs rhieey oor}! 


a a 
a Te. oo 4 wa Bouse, 


a | Gh Ched Hagetay) 'Ehat: 


Le . i¢ 4 ON 
- 7 ® 
: cian wae . 9g , ( 
7 


7 »>@ ws 
int : 
we . . 


bECRL Ty, ee | ea 


' 


_ 


THE UNIVERSITY OF ALBERTA 


APPLICATION OF THE FINITE ELEMENT PROGRAM 'ADINAT' TO ONE 
. DIMENSIONAL CONSOLIDATION 


FRANCIS PUI-KWONG YUNG 


A TABS S 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 
IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE 


OF MASTER OF SCIENCE 


DEPARTMENT OF CIVIL ENGINEERING 


EDMONTON, ALBERTA 


FALL 1984 


roN 
“r 


ae 


¢ 
# 


cS 


TO MY PARENTS 


j 

i ' im ! 

é : - 
i 
7 

‘ - iad 

i 

J 
5 » 
=. 


“py 
& i : 
‘6 
i iff ui ! ta 
j a . vm 
‘ ’ =u re } , na 
Ay : ar! chy te 
vase Le Ky Naas tae 
’ 1) ie ; Ay 
iW i 


ABSTRACT 


This thesis presents a study of the application of the 
finite element program ADINAT (Automatic Dynamic Incremental 
Nonlinear Analysis of Temperatures) to one dimensional 
consolidation. A few basic features of the program are 
reviewed, and linear and nonlinear consolidation analyses 
are discussed in detail. Specifically, a simple general 
governing equation is formulated for the primary use of 


predicting consolidation pore pressures. 


Applications of ADINAT in pore pressure predictions are 
illustrated through various sample analyses. Illustrations 
have included a comprehensive review of relevant researches 
to date. Within the scope of the present study, this 
Gissertation may serve as a state of the art manual for a 


generai treatment of one dimensional consolidation. 
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LIST OF SYMBOLS 


constant; Lagrangian coordinate, measure in the 
direction of gravity 


reduced coordinate, measure in the direction of 
gravity 


energy capacity; compressibility; reduced coordinate 
mode of variation of the modulus of elasticity 
Specific heat per unit volume 

distance from water table to soil layer 

void ratio; superscript indicating "element" 


Subscript for final value; bold type indicates load 
VeECron 


consolidation parameter 

load Ms 

total depth of compressible layer 
Subscript refers to nodal value or counter 


conductivity; permeability; bold type indicates 
VEOQEOr 


thermal conductivity 
Subscript denotes number of elements 
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moderotavariationt forec 
F 


coefficient of volume compressibility 


a constant; an outward normal to boundary surface; 
Subscript refers to values obtained with zero void 
ratio 


Subscript refers to the initial state 


pore water pressure; material constant 


energy flow rate per unit volume; fluid flow; 
applied loading; material constant 


rate of energy generated per unit volume 


rate of energy stored per unit volume 


boundary energy flow input to the surface S, 


radius 


Subscript for value at soil surface 


time 


EXCESS DOPCH DECSSUTEe 


volume flow velocity; subscript indicates vertical 
direction 


Subscript indicates water 


coordinate variable in ADINAT 


XV 


nh ‘ 
ah Dual a 


até 
iz iz 


rh, ai Wi: 7 


sna aati 
ages at Sh sans 

ae or qisbavod. es 
BNSF O20) Pete tices 


3 


ns tens bait 
oe O' 


Ej ihk, 3 


a4 


‘%o sitet 


anu lov oDAGy te 


re oa = 


SmipLow Sint! 14¢ pate Ye ovens to Bis 


4 
Lp . 
~ Hy pt - s - La at v 
-f- egsinjapeatror som: WELT NO IATS Yr SRVSG : 
pre | i 
= = ° . ; = . : ae 
au fDes 
: . 
e 
I7Gixzie. Stoz2. davebar oO? Fe essed | 
’ jae) 7] : , = 9) - 


betijtay S27 8orAanl 9ae ssedideas poten we BY em tov 
ra. ly es ianaatB 


io 


Co 


Ge 


Eulerian coordinate, measure in the direction of 
gravity 


a constant; cross-sectional area; bold type for 
interpolation matrix 


amplitude of sine function 


Terzaghi's excess pore pressure ratio; bold type for 


gradient interpolation matrix 

mass Or compressibility matrix 
material constant 

Slope of the e-log;,.0' line 
consolidation parameter, (1+e,/lte)?Cv 
material parameter, (dAe'+1)/(A+1) 
Slope of the e-log;,.k line 

material constant 

coefficient of consolidation 

Cv with constant m 
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given function; operator 
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G given function; operator 


H initial length of vertical drainage path; subscript 
refers to heat 


K Symmetric matrix; conductivity matrix 
K permeability variable in ADINAT 
KM refers to either the permeability or the 


compressibility 


L differential operator 
M dimensionless number 
M compressibility variable in ADINAT 
N Shape function; total number of solution time step 
ae] 
Q rate of energy generated 
Q' concentrated energy flow input 
oR. energy inflow rate 
out 
QO energy outflow rate 
o* rate of energy stored 
Q total nodal energy input vector 
Q energy generated at the system 
B 
Q energy stored at the system 
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energy transfer at the boundary surface 

externally applied nodal energy flows 

boundary surface of the problem domain; settlement 
boundary on which the unknown is prescribed 
boundary on which the unknown gradient is prescribed 


temperature; time factor; superscript denotes 
transpose of matrix 


time factor for bottom layer 
tame ofac tor, tC ot/h 
time «factor; (c'ky.it) /y 
W 
average degree of consolidation for the whole layer 


weight 


load ratio; space coordinate; time integration 
scheme 


See ratio; space coordinate 
unit weight 

average submerged unit weight of soil 
unit weight of soil solids 


variational operator 
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€ mode of variation of material properties; indicates 
"an element of' 


§ self-weight ratio 
i‘ average self-weight ratio 
n mode of variation of material properties 
n' material constant 
Gg unknown function 
0, prescribed unknown function 
6! unknown function at boundary point 
ae unknown function at the boundary surface 
6 vector of the unknown function 
r eigenvalue; material parameter, meo/(Co-meo) 
A! material parameter 
Li local consolidation ratio or relative degree of 
compression 
p rabior ofr Ce/Ck 
p Mater@aleeconstant,@1/Cci — 17Gk 
r 
E strain: convective coordinate 
p density 


x1x 


ere 


Hebe esisveqosq [sists mach (Ci 6 
< : — : 7 


jw. 


oa r 
: _ 
: ste 7eEnOS 
— 
is 
Te jn 2 
r Bet oWOoAT: 
ne 
stag Vise “2 35 #9 : ; 
p27 ue “Y Maes + 45°70 . 
oh Sd Le eH HOO zi 
{ - 7) ait vc 2 We rors - 
=e “Ais ~ £ 
Oo SeeD avcselieit 10 Of561 Tar7 
i] 
a 
Ph as ew J 
») 
a , ie anes ; ey 
a5 oO). gasetesoe Laz 
“ale ma ay 
» 
Je 


7 , ict Ts wi 
. agen 264 ay 
. . a neesye is 4) 


total stress 
effective stress 


initial effective stress in soil, before loading 
applied 


dimensionless time factor, Cvi-*t/H? 


Subscript for initial value at surface (before 
loading is applied) 


diffusivity 
Variable of logis (eo /ea ) 
15 
material parameter 
prefix indicates a finite increment 
Space coordinte denotes either @ or Zz 
a stationary functional 
Summation 
bounded problem domain 
characteristic number 
reduced thickness of soil layer 


Subscript or superscript refers to the ith layer 


Subscript refers to an arbitrary reference state 
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indicates the final value at layer bottom 


Superscript prime denotes coordinate differentiation 


quote mark refers to a ratio of the original value 


Superposed dot denotes time differentiation 


upper underscore indicates mean value or trial 
function 
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1. INTRODUCTION 


1.1 General 


The majority of phenomena arising in engineering 
and physics are often described by differential equations, 
especially by partial differential equations, and by 


boundary conditions imposed on the unknown functions. 


Seepage through porous media, heat conduction and 
torsion are three important phenomena governed essentially 
by differential equations of similar form, particularly the 
diffusion, Laplace and Poisson equations. These problems 
are often known as field problems, and the general governing 


equation is called a field governing equation. 


Some other examples of this type of problem include 
gravitational and irrotational fluid flow, surface waves on 
aaftdiuids natural frequency of footings, heat convection, 
electrostatics, magnetostatics and electrical currents. 


(Desai and Abel, 1972) 


ADINAT, a general purpose computer program which 
employs the parabolic differential equation of the heat 
conduction type, can be used to solve many of these field 


problems. 
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Since finding closed-form solutions for the more 
general cases can be difficult, various methods of computing 
approximate solutions have been developed. Among these, the 
energy method, or variational method, has gained 
considerable popularity. The finite element formulation, 
which is modified from the classical variational method, 


became widely used with the advent of computers. 


ADINAT employs the variational procedure in deriving 
the finite element equations. The high accuracy of the 
program and its incremental solution procedure are 
particularly applicable to problems with soils, where 
nonlinear constitutive relations are fundamental and 


Lnerinsic. 


Consolidation is the process of gradual reduction in 
volume of a soil mass of low permeability due to the 
dissipation of excess pore water pressure. The one 
dimensional consolidation condition occurs either with a 
large uniformly loaded area or with lateral physical 
constraints. One dimensional consolidation theory is 
extensively used in practice for pore pressure predictions, 


settlement and stability analyses. 


The conventional consolidation analysis is based on the 
classical Terzaghi theory. Despite the simplifying 


assumptions embodied in the theory, it is still widely 
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employed in practice due to its simplicity. It also permits 


extensions in a simple way to cover other situations. 


An example of an earlier extension is the 
nonhomogeneous problem. It is often approximated as 
discrete layers of unlike soils, in each of which the 
material properties are reasonably uniform. Continuous 
Spatial variation of material properties, is another simple 
method to simulate the varying material properties in soils. 
The latter is, however, a step closer to a nonlinear 


analysis. 


The term 'nonlinear' here refers to the condition of 
nonlinear constitutive relationships, or of material 


properties which basically vary with time. 


One dimensional nonlinear consolidation studies were 
developed by 1965. Barden and Berry (1965), Davis and 
Raymond (1965), Janbu (1965), Mikasa (1965), Raymond (1966 & 
1969), Gibson et a] (1967). Poskitte (1969) peDavis (1971) 
provided prominent nonlinear theories at that time. Gibson 
et al published a general theory removing many of the 
simplifying assumptions in the conventional analysis, and 
included the effects of both finite strain and self-weight. 
Extensions to this theory have been made by other 
researchers. (Poskitt, 1969; Simons and Beng, 1969; Monte 


and Krizek, 1976; DeSimone and Viggiani, 1976; Gibson et al, 
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Due to mathematical difficulties, recourse has been 
made by most nonlinear theories to further approximations in 


order to derive solutions. 


Finite element solutions are provided by ADINAT for 
various nonlinear analyses. Detailed discussions of the 
application of the program to nonlinear consolidation will 


be presented in this thesis. 


1.2 Scope of Research 


The basic aim of this thesis 1S a general study of 
ADINAT application in one dimensional consolidation. 


Necessary steps to achieve this goal are outlined below: 


1. Examination of basic features in ADINAT; 

2. Thorough review of consolidation theories; 

3. Assessment of available models applicable to ADINAT; 

4. Formulation of the appropriate model compatible with 
the illustrative theory. 

5. Simulation of ADINAT input and preparation POrmUTe 
results. 


6. Evaluation of the use of ADINAT in such analysis. 
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Steps 3 to 6 are repeatedly used throughout this thesis. 


The consolidation studies here are limited by the 
assumptions stated in the following chapters. The studies 
can be divided into two parts. The first part deals with 
linear and nonhomogenous problems, in which Terzaghi's 
Simplifications are generally applied. The second part is 
concerned with nonlinear analyses, where assumptions such as 
the validity of Darcy's law and constant load increments are 
accepted but stress history and secondary effects are 


neglected. 


Applications of ADINAT are illustrated through sample 
analyses selected from various areas of the program 
applicability. Predictions of pore pressures are the 
primary concern in these examples. The guidelines given 


below are followed in the development of the examples: 


a. Simplicity and brevity; 

b. Usefulness, in the sense that each example 
represents a certain type of formulation; 

c. Correctness in modelling, by presenting the required 
derivation and a comparison with published 
formulations and/or experimental results; 

d. Correctness in computation, by comparison with 


closed-form or published results where available; 
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e. Completeness; the main reasons for those not 
illustrated are: similarity, duplication, inadequate 
data, no practical solution exists or the analyses 


are outside the scope of the assumptions. 


A brief discussion of the governing equation and 
several features in ADINAT are given in Chapter 2. Chapter 
3 deals with the application of ADINAT to linear problems. 
Conventional analysis, as well as homogenous layered and 


continuous nonhomogeneous problems, are considered. 


In view of the variety of soil stratigraphy, loading 
and boundary conditions that exist in nature, uses of ADINAT 
in different physical situations will also be illustrated 


and given with the examples. 


A brief review of theories and discussions on the use 
of ADINAT in nonlinear analyses are presented in Chapter 4. 
A collection of equations are given for a general treatment 
of these analyses. These equations are tabulated for 
convenience in use and reference. Chapter 5 and 6 are 


devoted to examples of various published nonlinear theories. 
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2. BRIEF DESCRIPTION OF ADINAT 


The governing equation in ADINAT will be discussed 
below. Mass matrix or compressibility matrix analysis and 
the time integration scheme of the program will be briefly 
discussed. Stability and accuracy requirements in using the 


program are also given. 


2.1 Governing Field Equation 


The general field equation solved by ADINAT, for a 
problem domain idealized as an assemblage of finite 
elements, can be interpreted physically as a mathematical 
Statement of the energy equilibrium at the nodes of the 
system at any time, that is (Wilson, Bathe and Peterson, 


1974: Holman, 1976) 


OMe l set Q Oye 2e4;) 


in which 
Q'" +: rate of energy flow into the elements adjacent to the 


node and into the node from an external source, 


ed | Genera ois nastosang 
=< A hell > oa * | | i 


4 eel bistt painzsv0d 11s 


Tin i 
¢ 16% (TREIOA yd ier doe | ao léeuge ‘bisii Layered aT : 


Saints Bie se ah blast Aa’ ae (beist Saab! fhiisaod neta 
Pep damedtaiaat go 24 ett fest cetedty Bggemgreste: Sd miso e2nemess. 
any 35 sbi sia +g cot ae eigmesan ef7 36 snomeaede | 


noens2e¢ ding re) a sh agi?) emia yas Js meteya 
| = Tarver .Welok 4oTe) 


hu. 


c ae. pati aria ¢ a 


su ® 
ania we ON 


soaeanentay nen : 


<) sal 
act: 


g 
Q : rate of energy generated within the elements 
adjacent to the node, 
out 
Q : rate of energy flow out of the elements adjacent to 
the node, 
OF : rate of energy stored within the elements adjacent 


to the node. 


Let @ denote the unknown function or the state 
variable, which can be fluid head or potential, stream 
function, velocity potential, temperature, warping, stress 
function, electrical potential and so on. The energy flux, 
q-dz, 1s generally related to the gradient of @. In 
considering the one dimensional condition, a very general 
relationship in terms of the proportionality constant, 


termed conductivity k,, is: 


06 
q.dz Se ee hee C2225) 
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where v, is the flow velocity and p, is the density. This 
can also be expressed by similar equations as those shown 
below using the appropriate physical law for a given 


problem. For instance, from Darcy's law 
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Or rrom Fourier’ s law of heat conduction 


g2,dz = -k, — (2.4) 


where u is the fluid pressure, T is the temperature, k, is 
the hydraulic conductivity in Eq. 2.3 and the thermal 


conductivity in Eq. 2.4. Hence, Q'" can be determined by 
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q is the rate of energy generation per unit volume while q 
is the rate of energy stored per unit volume. Combining the 


above relations gives, 


—(k,—) =q -q C2559) 


Generalizing this equation in three dimensions yields: 
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where. ky pp ky, k, are the conductivities corresponding to the 


ey era rections « 


The associated boundary conditions are 


6 = 6, C2aaine) 


on part of the boundary S, of the domain, while 


kp — = Q’ eee) 


on part of the boundary S2 of the problem domain, and 


Cite so <= G6 Cees) 


n denotes the outward normal to the surface. 6,, q* are the 
prescribed values of the unknown and energy flow input to 
the boundary surface respectively. 8) is the initial value 


of the unknown 6. 
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2.2 Variational Method 


2.2.1 Variational Principle 


A variational principle specifies a scalar quantity 
expressed in an integral form (Desai and Abel, 1972; 


Zienkiewicz, 1977) 
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in which @ is the unknown variable or function, the 
Superscript prime denotes differentiation with respect to 
the carteSian coordinates, the number of primes is the order 
of differentiation. 2 is the bounded problem domain and § 
is its piecewise smooth boundary, F and G are specified 
operators. Il is called a functional, because it is a 


£UREtDONTOE other functions. 


The boundary conditions can be sorted into three 
classes: the essential, natural and mixed boundary 
Gonditions., For quadratic functionals, where the operators 
contain at most second order derivatives, essential boundary 
conditions are the prescribed values of the unknown function 
6. The natural boundary conditions correspond to the first 


order derivatives of the unknown function, 06/dn. n denotes 
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ONnewormtNescOOrdinate Girectionex, vy, Zz, or the direction tof 
the outward normal to the boundary surface. The mixed 


boundary condition is a combination of the above two. 
In the variational procedure, tentative solutions, or 


trial functions, 6, satisfying the essential boundary 


conditions are tried for a given problem, 


6~=6=C Nia, (2545) 


where N,; are base functions prescribed in terms of 
independent variables such as the coordinates x, y, Zz. 


Parameters a; are unknown constants. 


From all such possible solutions, the solution to the 
problem will be the one which makes the functional I 


Stationary, that 15s 
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6 = 6 + 86 (2.49) 


6 1S the variation notation, termed the 'variation in". 


From calculus of variations, 


él = — 6a, + — 6a> +... = — Ja C2 ate) 


Since the stationary condition of II in Eq. 2.16 is true for 


any arbitrary 6a, then 


pean (e27a9)) 


which yields a set of equations for the parameters a. Note 
that Wfethepfunctionalull isoquadratic ,qorether deravativessof 
the unknown function @ do not exceed second order, Eq. 2.19 
reduces to a set of linear equilibrium equations of the» form 


(Zienkiewicz, 1977) 
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Ka + £ = 0 (2.26) 


fF is the action vector, and the matrix K will always be 
Symmetric. Hence the parameters a, are easily computed, and 


an approximate solution for the problem can be found. 


The essence of the variational method is, therefore, to 
establish the functional II of the problem, choose a trial 


Sinereonss and find the stationary value #e¢t Il, that is 


This approximation process using trial functions and 
Stationary conditions is known as the Rayleigh-Ritz method. 
The finite element form of this method employs several 
Simple trial functions a, each defined only over an element 
of a discretized domain, instead of using a single trial 
function which spans the entire domain as in the classical 


method. This piecewise fit is the basic feature of finite 


element analysis. 


The appropriate functional for a given problem can be 
established in a number of ways. It can be derived 
vigorously using some theorems as guidelines, or obtained 
backward from the governing differential equations and 
natural boundary conditions with mathematical manipulation. 


Generally, there does not exist a unique functional for a 
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particular problem. A specific functional, however, can be 


applicable to a certain class of problems. (Bathe, 1982) 


2.2.2 General Formulation 


A linear system of differential equations can be 


written as 


Ge 
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where § is the vector of the unknown variable, and g is the 
known action function. If the differential operator L is 
linear and self-adjoint, it can be shown that the variation 


principle is given as (Zienkiewicz, 1977) 


T dt 
TT = f [(1/2)6 Lé+6 g]dQ + Boundary Terms (2222) 
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the superscript T stands for the transpose. 
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Applying this to Eq. 2.10 and using integration by 
parts, the functional of the field equation is then (Bathe, 


1982) 


Bac 
- § @(q -q )AN - Jf B8qtdsS - TeiQi (99-73) 
Q i 


> 


2 


where @*, @' are the unknown functions at the boundary 


surface. Q' is the concentrated energy flow input. 


The functionals in engineering problems usually possess 
physical meaning. (Zienkiewicz et a], 1966) The first 
integral of the above functional, may be interpreted as the 
rate at which energy is being dissipated over the entire 


domain, where,. 
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exact solution, 1S consistent with the principle of minimum 
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where A and B are the element @ and @-gradient interpolation 
matrices respectively. Assume the rate of energy storage is 


given by 


qt sneeé (ares) 


the proportional constant c is called energy capacity (or 
volume compressibility). A Superposed dot indicates time 
ditrerentiation.« sThe variation in’ @, 66, in°EQ. 2525 can be 
considered as a virtual quantity and the equation is 
analogous to the virtual work equation in structural stress 
analysis. Hence the system governing energy equilibrium 


equations can be written as (Bathe, 1982) 
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in which C is the finite element system capacity matrix, 
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m is the number of elements in the discretized problem 


domain. K is the conductivity matrix of the system, 


T 
K = = JS B kB aQ (223'3)) 
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Q is the total nodal point energy input, 
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Q fis  asvector*of externally* applied nodale energy flows. °C 
E 
and K are constant matrices in linear steady analysis. 


ou ee _ . 


m4 


Bie,” <a ai ay ‘ 
. : vi re i as a 
- ' 1) ane Gun Wea) 


Alig 
, Wee 


lg Ag ae ee 
am 7 ; 7 
yougai cetene Jaliq labo Legos ot ef ‘ 


os) 


i pt .o) 


2A 


For nonlinear or transient problems, C and K may be 


either @ or time dependent or both, as 


Cle OEO FER (Ont) nel Of td (2.38) 


where t denotes time. 


2.3 A Note on The Lumped and Consistent Matrix Analyses 


Two types of matrix analysis options, the lumped 
and consistent matrix method, are available in ADINAT. Some 
of their basic features are briefly described in this 


section. 


2.3.1 Consistent Matrix 


Application of the standard finite element method to 
the energy equilibrium equation produces the capacity matrix 
Ceinetherrorm of ‘Eq. 2.32. This matrix iis’ called consistent 
because it is derived from a procedure according to the 
variational formulation. (Desai, 1979) The resulting matrix 


is non-diagonal. The term consistent is thought to be 
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unnecesSary Since it is the natural outcome of the 


discretization process. (Zienkiewicz, 1977) 


The principle advantage of the consistent matrix method 
is to provide a more rational analysis. However, its use 
does not always lead to improved accuracy but always 
involves additional computational work compared to the 


lumping method. (Desai and Abel, 1972; Zienkiewicz, 1977) 


In the frequency analysis using the consistent matrix 
approach, more accurate mode shapes and upper-bound 


frequencies are obtained. (Desai and Abel, 1972) 


2.3.2 Matrix Lumping Approach 


An alternative way of establishing an approximate 
matrix C is by physical reasoning. As such, the matrix is 
assembled in an arbitrary manner by assuming that all the c 
values tributary to a node is lumped at that node, thus a 
diagonal matrix always results although no actual 
concentrated values exist. ThisS matrix is often referred to 


as lumped matrix. (Zienkiewicz, 1977) 


The resulting diagonal matrix in the conventional 
lumped matrix approach leads to a simple technique of 


formulation and solution. In addition, the lumped matrix 
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utilizes less storage space. This method is hence 
Significantly more convenient and economical for many 
computational processes compared to the consistent matrix 


approach. (Zienkiewicz, 1977) 


But if a coarse finite element mesh is used the result 
obtained may be very inaccurate. Surprising solutions may 
also result in some cases when high-order elements are 


employed. (Bathe, 1982) 


The lumping scheme, nevertheless, can often be 
constructed so that there will be no accuracy loss. 
Reseachers have found some matrix lumping techniques that 
improve accuracy. (Fried and Malkus, 1975; Hinton, Rock and 


Zienkiewicz, 1976; Zienkiewicz, 1977) 


The consistent and lumped matrix option can be found in 


ADINAT users manual input section II.2. 


2.4 Remark on the Time Integration Scheme 
2.4.1 General 


For transient analysis, a one parameter (a) family of 


one-step numerical time integration method is employed in 
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ADINAT. This method, as in the general incremental 
procedure, assumes that the unknown functions at time ty 
have been determined and that they are to be computed at 
Conerurtnt yawhere: t, is the time at time: n,;.nei0;, 5.6, 
N-1}, N is the total number of solution time steps and At is 


the time step interval. 


To evaluate the unknown functions at time t,+At, the 
energy equilibrium of the problem domain is considered at 
time t,taAt. From Bq. 2.38, (Bathe and Khoshgoftaar, 1979: 


Bathe, 1981) 


Gleaner OereK( 0.) tape, = OCe,) (218) 
De aG Avo O 7 Ae (2.40) 

6., = (1-a)6, + a6ns1 C22 7Ans) 

aera te (2.42) 


In the above equations t, denotes t,taAt, or (nt+a)-At. 
Subscripts n and nt+1 denote the time of t, and t,+At 
respectively. 9) is the prescribed value at the initial 


state. 
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Note that a iS a constant that ranges from 0 to 1. The 
value of a for a given problem is chosen to obtain optimum 


Stability and accuracy. 


Some known integration procedures such as Euler's 
foreward method, the trapezoidal rule, and Euler's backward 
method correspond to a values of zero, 1/2, and 1 
respectively. These options are shown in ADINAT users 


manual input section I1.5. 


2.4.2 Stability 


The properties of the integration method depend on 
the value of a that is used. They can be analysed by 
considering a typical single degree-of-freedom equilibrium 
equation, (Hughes, 1977; Bathe and Khoshgoftaar, 1979; 


Bathe, 1982) 


ES MG oe = i (2.43) 


where A can be considered as a ratio related to k andc. 
The characteristic equation of the recursion relation is 


obtained by using similar procedures from Eq. 2.39 to 
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Eq. 2.42 and solving for the unknown function @ at time 


Beta, 


Gn. oneS2e, 
TC t—@)yAEAy 
= —_—— |¢. (2.44) 
aAtN. 
where for stability, it is required 
[SP sett (2.45) 


From this, the following results can be observed: 

First, for 1/2<a<1, the above relationship is always 
satisfied for any At, which means that the integration 
method is unconditionally stable. 
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that is, the method is conditionally stable. 


Therefore, the conditionally stable methods, 
particularly the explicit scheme (a=0), will always result 
in instability as long as the time interval At exceeds that 
specified by the highest eigenvalue in the problem. The 
time saving computation in the explicit scheme may, however, 


compensate the requirement of using many small time steps. 


204.30AccUracy 


The accuracy of the Euler forward method (a=0), 
the trapezoidal rule (a=1/2), the Euler backward method 
(a=1) can be assessed using the Taylor series expansion. 
(Hughes,1977; Bathe and Khoshgoftaar, 1979; Snyder and 


Bathe, 1981) 


it can Se cinerds that, both of the Euler methods are 
first-order accurate whereas the trapezoidal rule is 
second-order accurate. The remaining a schemes can also be 
shown to possess first-order accuracy. Hence, the 
trapezoidal rule is the only second-order accurate method in 


the a-family of integration scheme. 
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The a values with the largest error would be a=0 and 
a=1. In the Euler backward scheme (a=1), the truncated 
terms in the time derivative equation obtained from Taylor's 
expansion are alternate in sign, which may result in 


Significantly better solution accuracy. 


2.4.4 Oscillation Effect 


The following important property should also be 
considered in the evaluation of the different integration 
schemes in the a-family method. Refer to Eq. 2.44 in the 
previous section, oscillatory solutions will result as the 
characteristic value § becomes negative (Zienkiewicz, 1977). 
Assume that [C°'][K] of Eq. 2.31 has a complete set of 
orthogonal eigenvectors with corresponding eigenvalues Aj, 
there will be no oscillation effect provided that (Wood and 


Lewis, 1975) 
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At <fAt =miy) max Ny (2.48) 
critical 


In each integration scheme for which ae<1, there exists 
a time interval beyond which § becomes negative or violates 
Eq. 2.47. Thus the Euler backward method (a=1) is the only 
one not showing this oscillation phenomena. The predicted 
solutions using this method converge monotonically to the 


exact solution as the time interval At decreases. 


Experience shows that, the time interval required for 
an accurate solution in the Euler forward method is 
considerably smaller than the value for stability. Even if 
the trapezoidal rule shows a higher order accuracy 
mathematically, it can give less accurate results than Euler 
backward method for larger time intervals. (Bathe and 


Khoshgoftaar, 1979) 


The critical time interval in Eq. 2.48 is related to 
the number of degrees-of-freedom. As the number of 
degrees-of-freedom increases due to the increases in the 
complexity of the problem or the base functions, the 
critical time interval becomes proportionately smaller. 
Therefore, it is impractical to eliminate oscillations by 


reducing time intervals. (Wood and Lewis, 1975) 
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The procedure with a=2/3, named the Galerkin scheme, 
has been shown to be practical as a compromise between 
accuracy and oscillation effects. Another optimal scheme, 
with a=0.878 ~ 7/8 was also suggested by Lambert (1973), 


using different reasoning. (Zienkiewicz, 1977) 


Wood and Lewis (1975) have illustrated that the 
Crank-Nicholson scheme with a simple averging process is 
effective in dealing with oscillations. They also observed 
that parabolic and cubic finite elements give better 
performance for accuracy and damping oscillations than 


linear elements, particularly with smaller time intervals. 


2.5 Discussion 


Eq. 2.10 is the basic governing equation in 
ADINAT. Various field problems can be analysed using the 
program by manipulating the variables k, 6, an and ae 
These may be consolidation, seepage, and problems in oil 


sands and cold region engineering, specifically in the area 


of geotechnical use. 


In the present work, illustrations will be given for 
B 
the case of one dimensional consolidation. q is then set 


equal zero and Eq. 2.30 is employed. There is no applied 
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mass flow loadings and hence they will be ignored in the 


program input. 


Notice that the stationary value of Eq. 2.23 not only 
includes the governing equation but also the natural 
boundary condition. Hence, for the ADINAT input, only the 
essential boundary conditions are required to be explicitly 
Specified, since the natural boundary conditions are 
automatically incorporated in the formulation when using the 
Variational approach. This is particularly effective in the 


analysis of the problems with complicated boundaries. 


Regarding the essential boundary conditions, nodes at 
impervious boundaries are deleted as degrees of freedom in 
the input stage. This simply means that the equilibrium 
equations Eq. 2.31 will not be considered for these nodes. 


(see ADINAT users manual input section IV) 


The neglecting of natural boundary conditions input and 
the mechanical generation of system equilibrium equations, 
as shown in Section 2.2, are the primary advantages of the 


variational approach as used by the program. 


The value of a used in the following examples is 1 
(implicit scheme), mainly for reasons of stability. Other 
choices of a values may provide better results, depending on 


the given problems. Examples showing the use of a values 
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other than 1 can be found in Part B of the program manual. 


The equilibrium iteration scheme used by ADINAT is the 
modified Newton-Raphson method. The number of time steps 
between equilibrium iterations, maximum number of iteration 
cycles and allowable tolerance for convergence are described 


in the program manual and input in control card I1.4. 
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3. LINEAR APPLICATION OF ADINAT IN ONE DIMENSIONAL 


CONSOLIDATION 


The application of ADINAT's linear capability in one 
dimensional consolidation will be illustrated in this 
chapter. Examples will include Terzaghi's classical theory, 


homogeneous layered and continuous nonhomogeneous problems. 


3.1 Terzaghi Classical Theory 


3.1.1 General 


The classical theory of one dimensional consolidation 
of clay deposits was developed by Karl Terzaghi in Use ite 


is founded on the following basic assumptions: 


1. The soil is homogeneous and is saturated. 

2. The soil grains and pore fluid are incompressible in 
comparison with the soil skeleton, and volume change 
is the result of changes in void ratio. 

3. The average fluid flow and compression are one 
dimensional. 

4. The surface applied load causes a uniform increment 
of vertical total stress throughout the soil layer. 

5. Darcy's law of permeability is valid at all 
hydraulic gradients and in the form of the average 
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fluid flow velocity related to the excess pore 
pressure gradient. 

6. The stress history of the soil stratum is ignored. 

(2) tne ySstructural viscosity effects, which control 
secondary compression or creep, are ignored. 

8. All initial total stresses are uniformly distributed 
throughout the soil stratum. 

9. The void ratio iS a unique linear function of the 
effective stress and is independent of time, that 
is, the coefficient of volume change is constant. 

10. The coefficient of permeability remains constant 
throughout the consolidation process. 

11. Strains, velocities and stress increments are small. 

12. The soil layer is comparatively thin such that the 
stresses arising from the self-weight of the solids 
and pore fluid are negligible compared with those 


applied. 


By equating the volume of pore water expulsion with the 
reduction in void volume of a soil element, Terzaghi 
established the following governing equation for one 


dimensional consolidation: 
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Oe) du 
(zit) Sos=(z ft) C31) 
Oz ot 


Cv 


where Cv is the coefficient of consolidation, u is the 
excess pore presSure, z 1s the vertical distance from the 
nearest drainage boundary to the element under 


consideration. 


Eq. 3.1 defines the immediate magnitude of the excess 
pore pressure u in an element at a specific distance z in 
the soil mass at time t, after an instantaneous application 
of a time-independent load Aq. Since the load causes a 
uniform stress increment, Ao , the initial excess pore 


V 
pressure is 


CzeOye= No “= AG (G22) 
Vv 


Hence the excess pore pressure u for a clay layer subjected 


to double drainage can be expressed as (Taylor, 1948) 
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where M=7(2m+1)/2 and m is any integer from 0 to ~. The 


boundary conditions are, for double drainage, 


UCOet I B=90 (3.4) 
wan et) £=a0 (ae5) 
for single drainage, Eq. 3.5 becomes 
du 
—(H,t) = 0 (3.6) 
Oz 
where H is the minimum drainage path. Note that Cv is 
generally expressed as, 
k 
Cv Cara7)) 
y m 
Ww Vv 


k is the coefficient of permeability and m is the 
Vv 


coefficient of volume compressibility, which is defined by 
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(3.55)) 


The term (i+e) in the above equation is assumed to be 


equal to (ite,) in this chapter. 


The average consolidation 


U for the entire stratum in linear theory is determined as 
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The above integrals are evaluated for the entire 


compressible layer. Substituting Eq. 


Eq. 3.9 gives 
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The analogies between the heat transfer equation in 


ADINAT and the consolidation equation can be seen below, 
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On oT 
Heat Transfer: k Se fite east) 
Bao 2 2 ot 
Om Ui du 
Consolidation: k =yn — (3% 12) 
OZ? wv ot 


k is the thermal conductivity. 
H 

The input and output parameters of ADINAT can be 
deduced from the above comparison. It is worthwhile to 
mention that in the program, the product of p and c forms a 
Single variable c", which is the specific heat per unit 


volume. 
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3.1.2 Sample Analysis 


Objective: 


To illustrate the application of ADINAT to Terzaghi's 
one dimensional consolidation theory and predict the pore 


pressure distribution in a simple soil layer. 


Problem Description: 


A simple single layer with drainage at both top and 
bottom is shown in Fig. 3.1. The applied loading and the 
material properties used in the analysis are given in the 


figure. 


Finite Element Model: 


For the finite element analysis, ten equal length one 
dimensional linear elements have been used. The consistent 
compressibility matrix option and the Euler backward method 


for time integration have been chosen for computation. 
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PERVIOUS 
CLAY LAYER 10 ecm 

PERVILOUS 

MATERIAL PROPERTIES: ONE DIMENSIONAL FINITE ELEMENT 
MODEL OF THE CLAY LAYER 
k = 1.0 cm/sec 
m, = 1.0 em-/9 NODE 
1 

Cv = 1.0 cm?/sec | eo Caen 
ows 1.0 g/cm “2 


AREA = 1 cm? 


Le 


Figure 3.1 Analysis of one dimensional consolidation with 


double drainage 
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Solution Results: 


Fig. 3.2 shows the excess pore pressure isochrones 
puedicted by ADINAT ati-time factors 0.1, 0.42and 0.9.9 Only 
one half of the layer is shown because of the symmetry 
condition. Closed-form solution of Terzaghi's analysis is 
also given in the figure and good agreement is shown. The 


ADINAT input datafile is enclosed in Appendix A.1. 


Remarks: 


The computer program used to compute the closed-form 


solutions is given in Appendix B.1. 


A listing of ADINAT input datafile for the case of the 
impervious base is included in Appendix A.1. The physical 
problem, loading and finite element model is the same as 
above, except the thickness now is 5 cm. The pore pressure 
results, as would be expected, are the same as the upper 


half layer of the previous case. 


User Hints: 


a. For the case of double drainage, the initial 
conditions are applied to all but the top and bottom 


nodes, at which the deleted (zero) pressure condition 
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code are assigned in the input. Only the pressure 
condition of the top node is ignored in singly drained 
layer. 

The output values of excess pore pressure can be 
obtained for some selected time factors Tv by 
deliberately manipulating the input values of the total 
number of time steps and time interval, which product is 
the solution time t. 

There are generally three approaches to obtain the 
required results. One of these is simply to solve for 
one Tv at a time. Several computer runs are then 
necessary for a range of required Tv values. This 
method will give a uniform degree of accuracy for all 
the results. 

The other two methods are to solve for all the 
required time factors in a single run, with either only 
one or various time intervals. The degree of accuracy 
for the solutions, obviously, will not be consistent for 


these methods. 
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3.2 Homogeneous Layered Problems 


Be 2. 1-General 


In the previous section, the consolidation of a 
homogeneous soil stratum is considered. Stratified soil 
Systems, however, are a common occurrence. Post glacial 
lake clays usually show laminated or layered structures, due 
to various sediments sizes during seasonal deposition. 
Problems of unlike soil layers are often encountered in 


earth dams, embankment buildings and landfills. 


In this section, the consolidation of layered soils is 
discussed and examples are given. Since the material 
properties are assumed constant in each layer, the problems 


are not considered as non/inear. 


A closed-form solution for two layers is given by Gray 
(1945) and for multiple layers by Barden and Younan (1969), 
Schiffman and Stein (1970). Numerical solutions are 
suggested by Abbott (1960), Luscher (1965), Christie (1966), 
Davis and Lee (1969). Experimental studies are provided by 


Raymond (1966) and Barden and Younan (1969). 


The governing equation, Eq 3.1, can be rewritten as 
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10.40 du 
Cv (256) = *—— (Z5't) rakes) 
OZ 4 ot 
eg) 
where the superscript 1 refers to the 1 layer. The 


initial and boundary conditions are the same as in the 
preceding section. In addition to excess pore pressures, 
the continuity of flow between adjacent clay layers requires 
that flow velocities be equal at the common layer 


interfaces. These conditions aresformulated as 


1 Ve 
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1 1 
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where Z is the distance from the surface to the interface 
1 

Separating thentthyand (1+1)** layer. The void ratio — 

effective stress mebatiaonship differs for each layer, as 
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Figure 3.3 Void ratio - effective stress relationship in 


layered soils 
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Examples, consisting of two, three and four contiguous 
soil layers with singly or doubly drained boundaries, are 
analysed below. Conditions of uniform initial excess pore 
pressure and time-independent loadings are considered here. 
Closed-form, experimental or numerical results will also be 


included to compare with the ADINAT solutions. 
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3.2.2 Two Layers 


Object ive: 


To verify ADINAT application to Gray's (1945) two layer 


consolidation problem. 


Brief: 


The first closed-form solution for a layered problem 
was published by Gray for two adjacent compressible strata. 
The impeded drainage condition was also considered in his 


paper. 


Physical Problem: 


Two adjoining layers of unlike compressible soils under 
a constant applied load are shown in Fig. 3.4. Boundary 


conditions and material properties are given in the figure. 


Finite Element Model: 


In the finite element model 20 equal length elements 
are applied for the clay stratum. Similar to the previous 


example, a consistent compressibility matrix option and 
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MATERIAL PROPERTIES: FINITE ELEMENT MODEL FOR 
THE TWO CLAY LAYERS 
LAYER I 
kng=C2.5R210°2 cm/sec NODE, 
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Figure 3.4 Consolidation analysis of two-layered problem 
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Euler backward method are chosen in the analysis. 


Results: 


The ADINAT results of excess pore pressure are plotted 
in Fig. 3.5 with Gray's solutions and close agreement is 
observed. The parameter Tb shown in the figure is the time 
factor for the bottom layer. ADINAT input data is given in 


Appendix A.2. 


Others: 


A Fortran program is used to compute the pore pressure 
according to Gray's formulation. The source code can be 


found in Appendix B.2. 


User Hints: 


Two approaches can be used for the data input in 
layered problems. The material properties for each layer 
can be input in ADINAT, either as different linear element 
groups in control card I, or as different sets of material 
properties in a single element group in control card X. The 
Same results will be obtained for these two methods. Both 
input methods are used in the sample analyses and they can 


be seen from the datafiles in Appendix A. 
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3.2.3 Multiple-Layer System 


Objective: 


To verify the ADINAT analysis of the Barden and Younan 
(1969) solution of one dimensional consolidation for a 


three-layered soil. 


Brief: 


Closed-form and experimental results are given by these 
authors for artifically prepared layered soils. The same 
problem here is also analysed by Desai and Johnson (1972) 


uSing a finite element scheme. 


General Descriptions: 


The singly drained three-layered soils and the applied 
loading are shown in Fig. 3.6. The material properties of 
each layer are given in the figure. Input data is in 


Appendix A.3. 


The finite element model is also included in the 
Fig. 3.6. In order to be comparable with Desai and 


Johnson's (1972) results, eight elements are used for each 
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MATERIAL PROPERTIES: 
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FINITE ELEMENT MODEL FOR THE CLAY LAYER 


LAYER TOP | mippteE | sorrom | 


NODE 1 fe) 17 25 
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Figure 3.6 One dimensional consolidation analysis of a 


three-layered clay system 
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layer. 


Results: 


Pore pressure disSipation at the interface between 
Derwent and Kaolinite is plotted in Fig. 3.7. The ADINAT, 
closed-form, finite element (Desai and Johnson, 1972) and 
experimental results are compared. Good correspondence can 
be seen among the first three solutions over the entire 
range of consolidation, and with the experimental results 


above a time factor of about 20. 
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Objective: 


To verify ADINAT predictions in one dimensional 
consolidation analysis of a four contiguous layer soil 
system with the Schiffman and Stein (1970) closed-form 


solution. 


Brief: 


These authors have studied a multiple layer system with 
a general set of boundary conditions as impervious, impeded 
and free draining. Unrestricted loading history is also 
considered and a closed-form solution is provided using 
Fourier mathematics. An example of a four layer 
overconsolidated glacial till deposit with doubly drained 
boundaries and constant load is given. The same example is 


illustrated here. 


Descriptions and Results: 


Fig. 3.8 shows the soil profile and the finite element 
model. The material properties given by Schiffman and Stein 
are tabulated in the figure. Input datafile is in 


Appendix A.4. 
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Figure 3.8 Consolidation analysis of a four-layered clay system 
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A set of excess pore pressure isochrones at 25%, 50% 
and 75% consolidation are computed using ADINAT. These 
results and the closed-form solutions from Schiffman and 


Stein are compared in Fig. 3.9, good agreement is observed. 
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Figure 3.9 Pore pressure isochrones of a four-layered system 


= 


60 


3.3 Continuous Spatial Variation of Material Properties 
3.3.1 General 


Another type of spatial variation, where the soil 
properties vary with depth in a continuous manner, is 
examined here. In this case, the strains are assumed to be 
so small that the permeability and compressibility are 
effectively independent of time but are continuous functions 
only of the space coordinate z. This type of problem is 
only a special case of the more general material nonlinear 
problem. The latter will be developed further in the next 


chapter. 


Values of permeability and compressibility at any depth 
vary greatly due to the diverse soil type, initial stress 
State, stress increment and past stress history. The range 
of possible variations in nature is so wide that only 
special cases of distributions which can be expressed in 
Special form are used in general analysis. Therefore, only 


these cases are discussed in the given examples. 


These known functions, however, can be useful to 
Simulate many real situations. A polynomial variation can 
be used to approximate the conditions in normally 
consolidated or overconsolidated clay layers with sufficient 


generality. For example, in the case of overconsolidated 
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clay, where great amount of overburden has been eroded away 
in the geologic past, the upper layer usually is more 
permeable and compressible than at depth. This is 


PePusSerated gin Fig’. sal0a. 


Fig. 3.10b shows the variation in a normally 
consolidated clay layer. The low permeability and 
compressibility at top is due to surface drying and the 
development of a hardened crust. As shown in Fig. 3.10c and 
Rig.) 3.100, ithe sinusoidal distribution is close to the 
condition of partly consolidated clay layers, in which the 
center part 1S more or less permeable or compressible than 
the edges. Examples with these distributions will be 


analysed below. 


The governing equation Eq. 3.1 of the conventional 
consolidation theory is extended to include spatial 
variations in the coefficient of permeability and the 
coefficient of volume compressibility as follows, (Schiffman 


and Gibson, 1964) 


(ee du du 
— —{k(z)—(z,t)} =m (z){—(z,t)} (3.16) 
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Figure 3.10 Polynomial and sinusoidal functions for material 


properties distribution 
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Section 3.1. 


Note that even if k(z) and m (z) vary in the same 
manner and thus make Cv(z) a sedeeene (see Eq. 3.7), 
Ege 3296 does not neducesto Eq: 3.1. This only takes place 
when k and m are separately constant. Eq. 3.16, however, 
is only an Ear ation to a general analysis of a thick 
layer with nonlinear properties. This will be discussed in 


the next chapter. 

Note further that values of average consolidation for 
nonhomogeneous problems are determined from the settlement 
ratios. The same applies to nonlinear problems in the 


following chapters. 


The polynomial relationship is expressed in the form of 


E 


KMUz) =e KMo 1272) ee) 


while the sinusoidal variations are given by 


KM(z) = KMg + A*sin(2z') (ent: 
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in which KM indicates either the permeability or 
compressibility, z' is the depth ratio z/h. The subscript ¢ 
denotes values at the top of the soil layer. The parameters 
e and n describe the mode of variation of k andm. A is 


Vv 
the amplitude of the sine function. 


Published analytical solutions can only be found for a 
limited number of particular cases. Edelmann (1953) 
considered the linear variation in the modulus of elasticity 
and Martins (1965) has analysed the case of a linear 
variation in permeability. Closed-form solutions are 
provided for both cases. Schiffman and Gibson (1964) have 
considered a series of cases for the spatial variations in 
permeability and compressibility. Known distributions such 
as polynomial and sinusoidal functions are used in their 


analysis. 


7 . 


He Bes he. . Ay, isla 
i ialiniedi eat jie 


; ee — 
are semonen ant aet 2 


- ft 
lia sein | ae ree alten 
ov : y 
ny . 
1 : a7 : 7 . ’ 7 : 
| | ) gnotawe due, 
ot bowel ed yite eso aROnmee 3 festt< fens bode dud. 
4 - L os Se i -*. 4 ait . 7 
; ‘wee af cal ; 
Es . nramwishie f : sfarstree =¢ tedmun simit 
ae ; if 


eae 
try) ce 


DAE: ve 4 Pati ) roma i) Doria f 2 ita He foe Sood 


rotted ae [epee ait 1pd-SBSS> I SLT 


ru ige mag’ - hae “si9 


shave at Ds He bas is2 lish fonue ‘a bee fsimoag 


“2 


rt 4 


1 Hy _ i 
in} 


65 


3.3.2 Variable Permeability with Constant Compressibility 


Objective: 


To verify ADINAT application to one dimensional 
consolidation with linear variation in permeability and 


contant compressibility. 


Brief: 


Martins (1965) studied the special case of linear 
permeability and constant compressibility. He assumed e=1 
in Eq. 3.17 for compressibility and use a polar coordinate 


transformation of Eq. 3.16, which becomes 


ara o-u du 
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Leet ns)! 2 Oe 2a1) 


kg 18 the value of the coefficient of permeability at the 
layer surface. This equation can be considered as the 
diffusion equation of an infinite cyclinder with radial flow 
and constant diffusivity, where Carslaw and Jaeger's (1959) 


solution is then applied. 


General Description: 


Fig. 3.11 shows the soil profile and the boundary 
conditions of the problem. The material properties and 
finite element model are given in the figure. Consistent 
compressibility matrix computation and Euler backward method 


are used. Datafile is included in Appendix A.5. 


RESUIitS: 


The pore pressure isochrone predicted by ADINAT at 50% 
consolidation is given in Fig. 3.12, and compared with 
Martins' closed-form solution. Good agreement is shown. 
The maximun discrepancy is about 4% at the bottom part of 


the curve. 
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Figure 3.11 One dimensional consolidation analysis with linear 


variation in permeability 
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Figure 3.12 Pore pressure isochrones: linear variation in 


permeability 
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Remarks: 


A Fortran program is used to compute Martins’ 


formulation. 


Appendix B.3. 


A listing of the source code is included in 
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Objective: 


To show ADINAT analysis of one dimensional 
consolidation with polynomial distribution in permeability 


and constant compressibility. 


General Description: 


Eq. 3.16 and 3.17 are applied in the analysis. The 
Parameters e€ and 7 estimated from Schiffman and Gibson 
(1964) are used. Fig. 3.13 shows the soil profile and the 
finite element model. Forty equally spaced one dimensional 
elements will be used for the clay layer. Euler backward 
method and lumped compressibility matrix are utilized in the 


computation. The ADINAT input data is given in Appendix A.6 


Results: 


The ADINAT output of excess pore pressure at 50% 
consolidation is compared with the solution of Schiffman and 
Gibson and presented in Fig. 3.14. Good correspondence is 


observed. 
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Figure 3.13 One dimensional consolidation analysis with 
polynomial variation in permeability 
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Objective: 


To show ADINAT application to one dimensional 
consolidation with constant compressibility and sinusoidal 


variation in permeability. 


General Description: 


Eq. 3.18 is employed with the amplitude A obtained from 
Fig. 13 in Schiffman and Gibson (1964) paper. The soil 
profile and material properties are given in Fig. 3.15. The 
finite element model and the time integration method are the 
same as the last example. ADINAT datafile is shown in 


Appendix A.7. 


Results: 


Fig. 3.16 shows the comparison of the resulting pore 
pressure isochrone at 50% consolidation and the solution by 


Schiffman and Gibson. Good correspondence is noted. 
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Figure 3.15 One dimensional consolidation analysis with 


sinusoidal variation in permeability 
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3.3.3 Variable Compressibility with Constant Permeability 


Objective: 


To verify ADINAT analysis of the variable 
compressibiltity and constant permeability model with 


EGelmann’s (1953) solution in one Gimensional consolidation. 


Brief: 


Edelmann assumed the linear variation of the modulus of 
elasticity, E, in a soil layer whence the coefficient of 


compressibility becomes a rational function of depth, 


Mecte— 9 1 / (ery 2 to (3322) 
Vv 


where og iS a small initial effective stress at surface and 
c' is a constant characterize the manner of variation of E. 
PUtting this into=bo. 3.16, Which cansbe cimplitired into a 


different form 
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Edelmann employed Fourier mathematics to obtain a 


Td 


Coa 230 


(3.24) 


(537,251) 


closed-form solution for the case of a single layer with 


only top drainage. 


Description of Problem: 


The singly drained soil profile and loading are 
presented in Fig. 3.17. The compressibility model and 


material properties are given in the figure. 


datafile is given in Appendix A.8. 
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Figure 3.17 One dimensional consolidation analysis with 


variable compressibility 
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Finite Element Model: 


The finite element model consists of 40 equally spaced 
elements. Consistent compressibility matrix option and 


Euler backward method is employed. 


Results: 


Edelmann's closed-form solution at 50% consolidation is 
compared with the ADINAT pore pressure results in Fig. 3.18 
and good agreement is shown. Maximum error of the pore 
pressure ratio occurs at the top part of the plot and is in 


ehemorden oc: 10F +; 


Remarks: 


A Fortran IV program is prepared for Edelmann's 
solution. The roots for the Bessel functions in the 
calculation can be found from Abramowitz and Stegun (1972). 
The required time for the analysis 1s computed by the 
program. A listing of the source code is included in 


Appendix B.4. 
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Objective: 


To demonstrate the ADINAT application in one 
dimensional consolidation with sinusoidal variation in 


compressibility and constant permeability. 


General Description: 


The soil profile and loading for the problem is 
described in Fig. 3.19. With the values e and 7 estimated 
from Schiffman and Gibson (1964), the compressibility model 
1s obtained and given in the figure. Other material 


properties are also included. 


The required time in the analysis is estimated by 
Simpson's rule in the numerical integration of the average 
settlement. A sample of ADINAT input data is given in 


Appendix A.9. 


The finite element model, time integration method and 
compressibility matrix option are the same as the previous 
example. These can also be found in the datafile in 


Appendix A.9. 
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PERVIOUS : ee yet 
a err 


CLAY LAYER 20 cm 
PERVIOUS 2 
My/My¢ 
COMPRESSIBILITY MODEL: 0 


My = Myo - A*Sin(mz"') 
Myo = 1.0 cm?/g 
eee aC eC. 7c] ' 


MATERIAL PROPERTIES: 


kK = 1.0 cm?/sec 


~ 
" 


a PO g/cm: 


FINITE ELEMENT MODEL FOR THE CLAY LAYER 


0.5 cm 


> - AREA = 1 cm? 
41 


40 EQUALLY SPACED ONE DIMENSIONAL ELEMENTS 


Figure 3.19 One dimensional consolidation analysis with 


Sinusoidal variation in compressibility 
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Results: 


The excess pore pressure isochrone at 50% consolidation 
predicted by ADINAT is compared with the solution by 
Schiffman and Gibson in Fig. 3.20, and close agreement is 


found. 
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Figure 3.20 Pore pressure isochrones: sinusoidal variation in 


compressibility 
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3.3.4 Geometric Variation of Material Properties with 


Constant Coefficient of Consolidation 


Objective: 


To illustrate the application of ADINAT in the analysis 
of one dimensional consolidation with variable permeability 
and compressibility but constant coefficient of 


consolidation. 


General Description: 


Fig. 3.21 gives the soil profile and the finite element 
model. The estimated values of the parameters for the 
material model are obtained from Schiffman and Gibson (1964) 
and given in the figure. The estimated time is found from 
Simpson's rule for the integration of the average 
settlement. The ADINAT data is included in Appendix A.10. 
Finite element model is the same as previous example and is 


also included in the figure. 


Results: 


The pore pressure predicted by ADINAT is presented in 
Fig. 3.22. This same problem is also analysed by Schiffman 


and Gibson (1964), close approximation is shown. 
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PERVIOUS 
CLAY LAYER 20 cm 
PERVIOUS 


k/ko, my/Myo 


{ 
MATERIAL MODEL: e 
€ 
k = ko (1 =o 7Z 2) 
€ 
My = Myo(1 + nz') 
z' 
kg = 1.0 cm/sec 
ive 1.0 em-/g 
No =18 
e =f1.05 1 
OTHER MATERIAL PROPERTIES: Yw Ce! et 6. g/em: 
Cv = 1.0 cm?/sec 


FINITE ELEMENT MODEL FOR THE CLAY LAYER 


0.5 cm 
+ - ea = 1).cm- 
NODE 1 41 


40 EQUALLY SPACED ONE DIMENSIONAL ELEMENTS 


Figure 3.21 One dimensional consolidation analysis with 


polynomial distribution in material properties 
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Figure 3.22 Pore pressure isochrones: variable permeability & 
compressibility with constant coefficient of 
consolidation 
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3.4 Discussion 


In most of the above examples, closed-form solutions 
are given for a qualitative evaluation of ADINAT 
predictions.sniThewhigh accuracy of the ADINAT results found 
in the examples in this chapter serves as a check on the 


reliability of the program when applied to such analyses. 


PiesaverageecOstsOfh Gach Tunais about thirty tower acy 
cents for the size of above problems. This indicates the 


economy in using ADINAT. 


Since the excess pore pressure u is the primary 
unknown, the finite element procedure guarantees the 
continuity of u, but not the pore pressure gradient, 3u/dz. 
The continuity requirement in layered problems, Eq. 3.15, is 
therefore not strictly satisfied for the case of linear 
finite element model in ADINAT. However, the accurate 
predictions in the above layered examples assure the 


applicability of ADINAT in these analyses. 


For layered problems with severe changes in material 
properties, however, higher-order finite element model may 
become necessary. A further discussion on this subject can 


be found in Desai and Johnson (1972). 
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For the analyses in this chapter, pore pressure results 
are obtained directly from the program output. Subroutines 
Or supplemental programs may be then required for the 
computation of other quantities such as average pressure 
dissipation or settlement. The settlement ratio becomes a 
more important measure for comparisons in nonlinear 


problems. 
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4, GENERAL STUDY OF ADINAT APPLICATION IN NONLINEAR 


CONSOLIDATION 


Nonlinear consolidation theories have been developed 
for decades and numerous studies can be found. The scope of 
discussions in this and subsequent chapters are in 
accordance with the first seven assumptions given in 


Securonrasi i t< 


Published nonlinear theories are reviewed first. The 
general governing equation is then derived, in order to 
facilitate the discussions that follows. It will be shown 


to embrace most nonlinear theories previously cited. 


From this equation, different formulations with the 
common major simplifications will be generated. They will 
be presented in a form suitable for use with ADINAT. These 
formulations form the basis for the numerical treatment of 
various nonlinear theories using the program. Applications 
of these ADINAT equations to existing theories will be 


illustrated in subsequent chapters. 
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4.1 Review of Nonlinear Consolidation Theories 


For many years, attempts have been made by researchers 
to modify the assumptions implicit in Terzaghi's theory to 
agree more closely with the physical conditions and material 
properties in real soils. The changes of permeability and 
compressibility in a mass of real soil are two significant 
departures from Terzaghi's classical theory, which result in 


a nonlinear theory. 


The first attempt to deal with these nonlinearities was 
by Schiffman and Gibson (1964), as illustrated in 
Section 3.3, where they considered the material properties 
as spatially dependent. Since the effective stress also 
varies with depth, the nonlinear behavior is then treated as 
a geometric nonhomogeneous problem. This, however, is only 
an approximate approach to the nonlinear stress - strain and 


Stress - permeability behavior. 


Permeability is known to vary with effective stress in 
a complicated manner. A widely used law is the linear 
relationship between the void ratio and logarithm of 
permeability. Raymond (1966 & 1969), Poskitt (1969), 
Poskitt and Birdsall (1971), Berry and Poskitt (1972) have 
used this relation in their theories. Other relationships 


are also proposed by Schmid (1957), Barden and Berry (1965), 
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Hwang (1966), Samarasinghe et a](1982), Juarez-Badillo 
(1983). Schiffman (1958) and Koppula (1970) are special 


cases of above. 


For the stress range normally encountered in practice, 
the usually accepted compressibility law is the linear 
relationship between void ratio and logarithm of vertical 
effective stress. Barden and Berry (1965), Davis and 
Raymond (1965), Mikasa (1965), Raymond (1966), Poskitt 
(1969), Wilson and Hwang (1969), all used this relation in 


developing their analyses. 


It is found in some soils that the coefficient of 
consolidation varies much less than either the permeability 
or the compressibility. Researchers such as Janbu (1965), 
Davis and Raymond (1965), Mikasa (1965), Davis (1971) assume 
this cofficient to be a constant in their formulations. The 
assumptions of constant coefficient of consolidation and the 
logarithmic effective stress - void ratio relation have also 
been shown from laboratory tests for normally consolidated 


clay, and where the secondary effect is ideally absent. 


Many of these theories, such as Barden and Berry 
(1965), Davis and Raymond (1965), Raymond (1966), have 
introduced nonlinear behaviour in thin layered soils 


although the restraint of small strain still persists. 
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Small strain analysis is based on the assumption that 
the deformation is small compared to the thickness of the 
consolidating layer. Therefore, it is not appropriate for 
an extremely compressible soil and only applicable to 
oedometer tests or soil samples sufficiently thin that 


self-weight effects can be neglected. 


The usual coordinate system employed in small strain 
theory is the Eulerian system where all measurements are 
referenced in space to a fixed datum. Settlement 
computation based on small strain theory, therefore, may 
probably be in error since the movement of boundaries is 


assumed in a fixed coordinate system. 


Finite strain analyses for thin layers have been 
developed by McNabb (1960), Mikasa (1965), Gibson et a/ 
(1967), Poskitt (1969), Wilson and Hwang (1969), Mesri and 
Rockhsar (1974), Lee and Sills (1979). The first presents 
the fundamental nonlinear equations without providing any 
solution, while the last uses a moving boundary approach to 
account for large strains. Mikasa assumed a constant 
coefficient of consolidation and used the compression strain 
as the primary unknown in his governing equation. 
Settlement analysis 1S given by Wilson and Hwang for 
different load ratios and permeability parameters. Mesri 
and Rockhsar have carried out an extensive study of 


nonlinear behavior, but failed to account for the influence 
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of soil depth. 


The general theory formulated by Gibson et a] (1967) 
reduced many limitations in the conventional analysis. It 
represents a basis for various future researchers. Studies 
using thin layer approximation have been made by these 
authors and other researchers. Their original equation with 
logarithmic material laws and further assumptions is solved 
analytically by Poskitt (1969) using a perturbation method. 
A closed-form solution is found by Simons and Beng (1969) 
with other approximations. Comparison with experimental 


results is presented by Burland and Roscoe (1969). 


Experimental investigation and numerical solution of 
Gibson et al's (1967) general formulation is suggested by 
Monte and Krizek (1976). In their>analysis, the liquid 
limit state is taken as the datumn for strain measurement. 
Linear relationship between constrained modulus and strain 
is employed, and the permeability is expressed as a second 


degree function of void ratio. 


In engineering projects such as large embankments, 
reclamation filling and groundwater lowering, very deep and 
soft clay layers may be involved. Large settlements may 


occur and effects of self-weight cannot be ignored. 
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Raymond (1969) and Davis (1971) have considered effects 
of layer depth in the context of small strain. Numerical 
solutions for Davis (1971) is suggested by Viggiani (1973). 
Variable load increment is also considered by Raymond (1969) 


and Viggiani (1973). 


Studies accounting for both the finite strain and 
self-weight effect are given by Janbu (1965) and Gibson 
et a] (1981). The former gives a general formulation 
Similar to Gibson et a] (1967), and particular cases are 
considered. The latter uses the same general formulation as 
in an earlier paper (Gibson et a], 1967). Numerical 
solution for Gibson et a] (1981) is provided by DeSimone and 
Viggiani (1976) using assumptions of constant coefficient of 


consolidation and logarithmic compressibility law. 


Despite the fact that Gibson et al's (1967, 1981) 
governing equation allows the general consideration of many 
factors influencing consolidation, other approximations are 
maeOEDOrated in order to arrive at a solution. This is a 
serious drawback to most extensive nonlinear theories, where 
mathematical difficulties usually arise for a general 


Solution. 


The advantages of ADINAT in such a Situation is 
apparent since its finite element procedure, accompanied 


with the nonlinear solution strategies, provide a powerful 
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tool for solving any form of expression for material laws. 
Since different soils are subject to different forms of 
material laws, ADINAT allows the choice of the most proper 
constitutive relations based on experimental evidence rather 
than being restricted to any particular form due to 


analytical convenience. 


4.2 Derivation of the General Governing Equation 


The one dimensional consolidation of a simple clay 
layer with drainage at both top and bottom is used for the 
illustration here. (Gibson et a], 1967) The soil profile 
and coordinate systems employed are given in Fig. 4.1. 

Fig. 4.1(a) and (b) show the Lagrangian system and refers 
all events to the initial (zero time) state. The convective 
coordinate, on the other hand, always moves with the actual 
material configuration. This is described in Fig. 4.1(c) 


and.¢d). 


TO further earity the following derivation, a 
pictorial description of the stress distribution in the soil 
layer in the convective coordinate system is shown in 
Fidimése2. } In the, following derivation, the first seven 


assumptions stated in Section 3.1 are still valid. 
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Figure 4.1 Coordinate system: (a),(b) Lagrangian: initial 
state, (c),(d) Convective: consolidating state 
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Referring to the convective system ~ in Fig. 4.2, the 


total stress o can be expressed as 


CUE) -oyerdeo(tleatto, + f 7, (£)dE + Ag on 
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dis the distance between the water table and the soil 
layer, S(t) is the settlement at time t, y,; and y are the 
unit weight of soil and water respectively. nie ferenbiaean 
Portuolly withsrespect to@t, £ = Hlayt), and trantorming 


into the Lagrangian system by the chain rule, a vertical 


equilibrium equation can be obtained 


ne ee ee (4.2) 
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determined, 


Le ae (4.4) 


The net inflow rate of fluid weight is 


3 eyw 
q -(q tdg ) = - —[(v -v,)——]-d&dadB (4.5) 
E E E o£ W i+e 


and the rate of change of fluid weight can be given as 


-dédadB] (4.6) 


By continuity, equating the above two equations and putting 
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W 0 eyw O& 
] = - —[——-—] (4.8) 
da Ww 1+e dt ite da 


By Darcy's law, the fluid flux is proportional to the 


excess pressure gradient as 


aay a? ee (4.9) 


Substituting Eq. 4.4 and Eq. 4.8 into the above 


equation, one obtains 
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0a W da O& dt ite da 


From Eq. 4.2, 4.7, 4.10 and the effective stress 


principle, 
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Pag Oy 8) a ) 


gives 


W 
[(iteg)— - (y.-y )J} = - _— (4202) 


Note that eo iS not a constant but varies throughout the 
Soil depth. Eq. 4.12 is a general governing equation with 


effective stress and void ratio as the primary unknowns. 


4.3 Overview of ADINAT Application 
4.3.1 General 


One dimensional nonlinear theories can be classified 
into four main categories: small strain - thin layer, small 
Sthainise¢thickslayer; finite strain - thin layer and finite 
Strain - thick layer analyses. Applications of ADINAT to 


each category will be illustrated. 


Mathematically, small strain requires the volume at any 


time t, (1+e), is equal to the initial volume (ite,). Thin 
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layer analysis does not allow for the total stress variation 
with depth due to the self-weight of the soil mass, which 
means y; = y . This first and second limitation 

W 


individually is not imposed on the finite strain and thick 


layer analysis respectively. 


General formulations used in ADINAT for each of these 
categories are derived below. Their uses will be discussed 


latter. Examples are provided in the following chapters. 


4.3.2 Formulations of ADINAT Equations 


(a) Small Strain and Thin Layer Analysis 


As discussed previously, (1lt+te) equals (1+e,) and y, 
equals y are presumed for small strain and thin layer 
Br aiee Further, the Lagrangian coordinate @ can be 
replaced by the Eulerian coordinate Z, since the initial 


configuration holds for all time. Using these 


Simplifications in Eq 4.12, 
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which can be rewritten as 
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(4.14) 


(4,15) 


(4.16) 
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Utapt) =foo4+ Ade- of (z,e) (4.18) 


Differentiating partially with respect to z and t gives 


du d0° 
— = - — (40319) 
OZ OZ 

and 
du 00° 
SS (4.20) 
ot ot 


Substituting these equations in Eq. 4.15 and Eq. 4.17, 
together with the conventional definitions of Cv, Eq. 3.7, 
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fe) de de 
—|o—]-— (4.22) 


(b) Small Strain and Thick Layer Analysis 


Assume small strain, (ite)=(1+e,.) and thus replace @ by 
the Eulerian coordinate Z in Eq. 4.12, which can be 


rewritten as 


—{k— - Gigs pt eee — P5223} 


In thick layer analysis using the Eulerian coordinate, 
self-weight of soil mass is to be accounted for in the 


effective stress relationship as 
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Phe Subscript ¢ refers to the initial value at soil surface. 


Differentiating partially with respect to z and t give 
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== Ss —_- (4.26) 
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Substituting Eq. 4.14, 4.26 and 4.27 into Eq. 4.23, gives 
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(c) Finite Strain and) Thin Layer Analysis 


For finite strain and thin layer formulation, setting 


W¥eresy jiinsEo..4.N2f obtains 
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3 1+@5 Oa° yw de 
Bee liek poll) YE Se aes (4.29) 
da 1+te da 1+@, Ot 


Assuming Eq. 4.14 and using the definitions of Eq. 3.7 and 


Eq. 3.8, Eq. 4.29 can be expressed as 
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in thin layer analysis. oo refers to the initial effective 
Stress in soil. Differentiating partially with respect to a 
enatGeastitplar tovsection (a), and substituting into 


Eq. 4.30, gives: 
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— J=y — — CaS 5.) 
0a 1+e da w 1+@o ot 


(d) Finite Strain and Thick Layer Analysis 


To account for both large deformation and self-weight 
in finite strain and thick layer formulation, the convective 
system should be considered. The effective stress relation 


is 
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Substituting Eq. 4.36 into Eq. 4 
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(40535) 


(4.36) 


(4.37) 


(4.38) 


into Eq. 4.38 and rearranging 


f) 1+e, du mv(1te) du 
—I[ k — ] = y ——— — (4.39) 
da 1+e da Ww 1+€o5 ot 


Since (1+e,.) varies with depth, another reduced coordinate Db 
may be used. This coordinate may be physically interpreted 
as the volume of soil solids lying between the datumn and 


the point being considered, that is (McNabb, 1960) 
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The mathematics of the effective stress relation now becomes 


a little simplier, 
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Stipstituce Bo. 4.41) into Eq. 4.39 produces: 


—[{ — — ] = ym (1te)— (4.43) 
eb pire tab wv ot 


In reduced coordinates, the soil layer thickness is reduced 
to £, which is obtained from the integral, Eq. 4.40, 


evaluated throughout the depth. 


4.3.3 Remark on the Use of ADINAT Equations 


A summary table of the above ADINAT equations is given 
in Fig. 4.3. The equations are renumbered from A? to AQ in 
the figure for the convenience of later reference. 

Equations A? and A2 are Eq. 4.21 and Eq. 4.22 respectively. 
A3 1S an obvious outcome from Eq. 4.21 when the parameter m 
1s assumed constant. AG refers to Eq. 4.28 while A5 is a is 


direct result by assuming m constant. A6, A7, A8, AQ are 


corresponding to Eq. 4.33, 4.31, 4.39, 4.43 respectively. 


Notice that even if the forms of equation A7 and Ad are 
Similar, their uses are different. For the application of 


ADINAT, A4 must be evaluated as described in Chapter 6. 


The purpose of Fig. 4.3 is to facilitate the users 


Study of nonlinear consolidation using ADINAT. The 
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equations given in each category, although in different 
forms, are equivalent. This will further assist the use of 
the table. The table can be used in many ways. Examples 


will be shown in the succeeding chapters. 


In these equations, Cv and m are just notations 
presented in the conventional ths They are different from 
the mean values Cv, m obtained from the conventional tests, 
and can be Biber ituved by any appropriate expressions 
accomudang "oO. their jibasic definitions, Eq. s3s7 and) 3.8. 
Unless otherwise stated, the material parameters can be 


CONStancs Ory Variables. 


Since the logarithmic material relations are so widely 
accepted that a table of the finite strain - thin layer 
analysis with this approximation is given in Fig. 4.4. 
Several letter symbols adopted in the figure are defined 


below. a iS the commonly used load ratio, 


(aes One) 0 n (4.44) 


and B is the permeability ratio, 
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B =k a Ko (4.45) 
i 


This ratio is always less than 1 in monotonic compression as 
in the illustrated examples. The constant v 1S a ratio of 


the material coefficients 


Dea COM ECK (4.46) 


wierencc ris the slope of theJe=loq;oo line, and Ck vissthe 


Svoperotethe e-log,,k line. 4 *isedefined in a Similar way: 


po=—- (4.47) 


uw is termed local consolidation ratio and with a value 


Papgese trom 0. tom. ft 1s given by 


uo= —— (4.48) 


The application of this variable simplifies the equations as 
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well as their use, as may be seen from Fig. 4.4. If the 


logarithmic effective stress - void ratio relation is 


assumed, the excess pore pressure ratio can be determined 


from was follows 


a” (4,49) 


Thestotal tinal strain, ~-§ , 1S expressed in the form of 


eee (4.50) 


Cvi is the initial consolidation ratio defined as 


Cvi koa, Gives) in 10ms y=Ce 


W 


keoeo lina / yt (4.751) 
Mie he 


Another approximation utilizing a more general material 


law used by Janbu (1963) and Koppula (1970) is 
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de de op 
voane | — | _ (4,52) 
do’ do 06 
and 
k k OG <q 
—={[—] (4.53) 
1+e ite 4, a6 


These are also given in Fig. 4.4. The constants p and q 
range between -1 and 0, the subscript 1 denotes a reference 
State. The logarithmic constitutive relations just 


mentioned are special cases of these equations. 


Since average values of material parameters are often 
obtained from conventional laboratory tests, the 
formulations uSing these parameters are also included in 
Fig. 4.4. The values of k, m and Cv are assumed to be 
obtained independently from a hes Note that Cv may not be 
equal to k/(y -m ). 

woo 

For the same reason as above, sets of different forms 
but mathematically equivalent equations are given in 
Fig. 4.4. Derivations of these equation are long but 


Similar to those in Fig. 4.3 and will not be shown. 


Equations in the figure are numbered from B1 to B12 for 
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future reference. 


All of these ADINAT equations in Fig. 4.3 and Fig. 4.4 
are presented in a clear manner. for ADINAT input, and are 


given in the following form: 
) 06 a6 
— | K — | = M— (4,54) 


where x, 6, K, M represents the coordinate, unknown, 
permeability, compressibiltiy variable respectively. To 
illustrate, whatever expression is entered in the position 
'K' will be input as the permeability variable in the 


program. The same applies to the rest of the parameters. 


Equations with an unknown variable o° are primary for 
obtaining pore pressure results, while those with an unknown 
e are mainly for computing deformations. Both of these, 
obviously, can be used for either purpose with supplemental 


programs to reduce the original output. 
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4.4 Discussion 


It can be seen from Fig. 4.3 that pore pressure 
can still be maintained as the primary unknown. Hence the 
Strain or void ratio is no longer a privileged variable in 
nonlinear analysis, as is believed by some researchers. 


(Mikasa, 1965; Janbu, 1965; Gibson et a], 1967 & 1981) 


The difficulty with the general governing equation, 
Eq. 4.39, 1S that the settlement cannot be assessed 
directly. On the other hand, Gibson et a/'s (1981) equation 
may be useful for computing soil settlement, if a solution 


can be easily obtained. 


Regarding the use of the general equation, Eq. 4.39, 
inconvenience may seem to arise in the computation of soil 
thickness £, or the evaluation of the integral Eq. 4.40. It 
is believed that this cannot be avoided in any analyses 
dealing with both finite strain and self-weight effects. 
These two quantities, however, are readily and accurately 


computed by simple computer-based techniques. 


The general governing equation A8 or AQ (Fig. 4.3) is 
more useful than Eq. 4.12 since solutions are readily 
obtained by ADINAT for the former, while no general solution 
is yet available for the latter. Therefore, A8 or AQ 


represents the working governing equation for a general 
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treatment of one dimensional consolidation. Note again that 
both of these equations are basically equivalent to 


Bqtud?i2. 


As mentioned before, most theories can be recovered 
from the proposed ADINAT formulations. Examples are given 


below. 


Terzaghi's classical theory is an obvious outcome of 
ADINAT A1 when k and m are assumed constant. The governing 
equation found by ee and Raymond (1965) is equivalent to 
A2 with constant Cv. Barden and Berry (1965), Raymond 
(1966) are particular cases of Al where logarithmic material 


relations are used. 


Analyses of soils with large depth by Schiffman and 
Gibson (1964), Raymond (1969), Davis (1971), Viggiani (1973) 
are special cases of ADINAT small strain - thick layer 
formulation, AG. In comparison with A8, Schiffman and 
Gibson's (1964) general equation is still deficient in the 


Contexteot thick layer formulation. 


McNabb (1960) has given an equation similar to A7. The 
equation given by Poskitt (1969) can be obtained from A6 
with the logarithmic material properties and change of 
variables. It is identical to equation B70 in Fig. 4.4. 


Gibson et a] (1967), Simons and Beng (1969), Wilson and 
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Hwang (1969) are particular cases of AZ. Notice that the 
parameter Cv mentioned in the paper by Gibson et a/ is 


actually a mean value Cv in which (1+e)=(1+e,) is implied. 


me general governing equation Eq: 13 in Gibson ef a] 
(1967) or Eq. 19 in Gibson et a] (1981) is basically 
identical to A8 or AQ respectively. Solutions are provided 
for A8 and AQ using ADINAT but no general solution has yet 
found for Gibson et a/'s equation. Numerical solutions are 
given, however, only for an approximate form of the latter. 


(DeSimone and Viggiani, 1976; Gibson et al, 1981) 
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5. APPLICATION OF ADINAT IN NONLINEAR CONSOLIDATION OF THIN 


SAMPLE 


5.1 General 


In order to verify the proposed equations in Fig. 4.3 
and show their applications, a series of examples, most of 
which are taken from published papers, are presented. Thin 
layer formulations, which most of these papers considered, 
arepdiscussedsin-this*chapter. “Applications of! thick® layer 


theories will be given in the next chapter. 


Another purpose of these examples is to show the 
correlation among the mathematical equivalent counterparts 
in the same category. Therefore, as many as possible will 
be used in the same example. Applications of the equations 


in Fig. 4.4 will also be illustrated as well. 


ADINAT equations in Fig. 4.3 and Fig. 4.4 are referred 
EOmasw/s 66... candspias.. in the followind é@xamples..eiteas 
worthwile to note again that equations in Fig. 4.3 are 
general equations. That is, they can be used for any 
material laws. Equations in Fig. 4.4 are special cases of 
Bieumnrtes strains thin layer formulations in Big. 4.3, 
with the assumed constitutive relations shown at the top of 


each sub-category. 
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In these two figures, equations in each set are 
mathematically equivalent and are presented in a simplified 
manner. The choice among them is a matter of preference or 


may depend on the primary output of the unknown variable. 


For all examples illustrated here, unless otherwise 
specified, one dimensional linear finite element model is 
used, twenty equally spaced elements are assumed for the 
compressible layer. Ten solution time steps are assumed. 
Lumped compressibility matrix option and the Euler backward 


integration method are employed. 
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5.2 Examples of Samarasinghe-Huang-Drnevich Study 


Object ive: 


To show ADINAT application in Samarasinghe et a] (1982) 


=e 
tudies 
—— SS . 


ADINAT Nonlinear Equation: At, A6, B4, BS, B6 


These authors have found that the permeability relation 


2 
k = Ck —— (S21) 
1+e 
where Ck is either a constant or it may vary with void 
ratio. They also assumed a common linear void ratio - 
logarithm effective stress relationship. 
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Description and Analysis: 


An example of greyish sandy clay with the exponent 
n=5.2 and Ck=6.58x10°*m/day obtained from Samarasinghe et al 
ms iliustrated here, Fig. 5.1 shows the soid profile, 
loading and boundary Berea on ee Material properties and 


the finite element model are also given in the figure. 


Both small and finite strain analyses are performed for 
comparison. Since a logarithmic compressibility 
relationship is included, equations in Fig. 4.4 can be 


employed. Datafiles are given in Appendix A.11. 


Results: 


ADINAT prediction of pore pressure distribution at 50% 
consolidation is plotted in Fig. 5.2. Solutionsiwith A6 and 
B4 are coincident since they are virtually the same 
equation. The small discrepancy among the results of B4, B5 
and B6 may be due to the difference in numerical input into 
the program. The finite strain analysis in this example, is 


only little different from the small strain analysis. 
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Aq = 47.9 kN/m? 
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PERMEABILITY MODEL: COMPRESSIBILITY MODEL: 
e" Cra loguos 
ka=i Ck 
tae Gow, Ora16 
Ck = 6.58x10-‘ m/day 
i. =. 2 
OTHER PARAMETERS: Co = 0.622 
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= 9.81 kN/m:; 
se / 
FINITE ELEMENT MODEL FOR THE CLAY LAYER 
0.0985 cm 
7 - AREA = 1 m? 
_— to 
NODE l 24 


20 EQUALLY SPACED ONE DIMENSIONAL ELEMENTS 


Figure 5,1 One dimensional consolidation with nonlinear 


material properties: Samarasinghe et al, 1982 
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Figure 5.2 Pore pressure isochrones: thin layer analysis 
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5.3 Barden and Berry Analysis 


Objective: 


To show the use of ADINAT in a Barden and Berry (1965) 


analysis. 


ADINAT Nonlinear Equation: (A) Al 
(By Al ob i be moo neeO 


Brief: 


These authors have given a small strain theory 
utilizing logarithmic material laws. They further 


approximated the permeability relationship as 


kuetkalidn’ um) C5a2) 
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in which k is the value of k when u=0. 7' and n are 


constants. Their equation is given as 
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They resorted to finite difference method to solve the 


nonlinear equation generated. 


Description and Analysis: 


(A) An illustration of Barden and Berry's formulation 
with an assumed soil profile and material properties as 
given in Fig. 5.3 is presented. Values of 1 and 1/2, as 
referred to by Barden and Berry, are used for the value of n 


DMPC cs 


(B) Results of a laboratory test given by Barden and 
Berry iS analysed by ADINAT. Both small strain and finite 
Strain analyses are illustrated.’ Bi, B2, B3 and BG in the 
finite strain formulations are included for comparison. The 


loading and material properties are shown in Fig. 5.5. 
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Figure 5.3 One dimensional consolidation with Barden 


(1965) formulation 
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Figure 5.4 Pore pressure results in small strain - thin 


layer analysis 
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Results 


(A)Pore pressures predicted by ADINAT are plotted in 
Fig. 5.4. Not much difference can be seen between the cases 
of n equals 1 and 1/2. This has also been observed by 


Barden and Berry as well. 


(B)The pore pressure results given by ADINAT at 50% 
eonsolidatron are plotted in Fig. 5.6. “The finite strain 


solutions are found close to the small strain analysis. 


The discrepancies among the finite strain results are 
due to the fact that different parameters are used. Note 
that the numerical value of Cv used in equations B2 and B3 
is not exactly equal to k/(y -m ), whereas the parameters k 
and m are used in B7 and BG. 

Vv 

Fig. 5.7 compares ADINAT predictions, Barden and 
Berry's solutions and the test results. The experimental 
settlement curve is fitted to the ADINAT curve to obtain the 


required time, which in turn is used to determine the pore 


pressures. 
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20 EQUALLY SPACED ONE DIMENSIONAL ELEMENTS 


Figure 5.5 Nonlinear consolidation analysis with experimental 


results 
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Figure 5.6 Pore pressure distribution in thin layer analysis 
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Figure 5.7. Pore pressure dissipation: smal] strain - thin 


layer analysis 
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5.4 Davis-Raymond Theory 


Objective: 


To verify ADINAT small strain - thin layer formulation 


with Davis and Raymond theory (1965). 


ADINAT Nonlinear Equation: Al, A2, B2, B3 


Brief: 


These authors have given a small strain theory for thin 


layer and the governing equation 1S expressed as 


On iy 820 1 du Iy coe 
-Cv[ — + (— )?(—)? J =- — (5.4) 
OC nds 0 OZ ie Gye 
It can be rearranged into a different form 
S 80éw OW 
Cv = — (5 25)) 
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with 


@) FS logio an C5r os) 


The general solution can be determined through 


B 
Cire 0n l (O57 oe CSive) 
f f 


B is Terzaghi's excess pore pressure ratio. Eq. 5.4 or 5.5 
1S a special case of equation A2, and Davis and Raymond have 
assumed the logarithmic compressibility law and that the 


Ceoefticient of consolidation is constant. 


Description and Analysis: 


The soil profile, loading and material properties are 
Given in Fig. 5.8. Finite strain analyses are also carried 
out with equations B2 and B3. Eq. 21 in Davis and Raymond 
(Eq. 5.5 in above) is directly solved using ADINAT as well. 
Terzaghi's theory is also analysed by the program for 


comparisons. Datafiles are attached to Appendix A.13. 
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Figure 5.8 One dimensional nonlinear analysis with constant 


coefficient of consolidation 
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Results: 


ADINAT) pore pressure results at Tw=0.4 and; with a load 
ratio of 4 are presented in Fig. 5.9. The small strain 
formulations agree well with the closed-form solutions using 
Eq. 5.7. Solutions with A2 and with Davis and Raymond's 
Eq. 2Zivare almost identical» as they should.be... The finite 
strain solution shows a lower prediction while Terzaghi's 


theory gives an even lower result. 
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Figure 5.9 Pore pressure isochrones: thin layer analysis 
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5.5 Raymond Theory 


Objective: 
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To illustrate ADINAT application to Raymond's (1966) 


theory. 


ADINAT Nonlinear Equation: Al, A6, B7, B8, BI 


Brief: 


Raymond has developed a small strain theory which 


includes both variations in permeability and in 


compressibility. His general equation is given as 


in pe avesconstants 
Ene Soil surface, or 
equation is the same 


logarithmic material 


du aru 02434) Comaa® 
cs ] = ———— — (5.8) 
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as equation A? with the assumed 


laws. 
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Description and Results: 


Fig weoeOmeshows the soil, profile, and material 
properties. The loading and finite element model is given 


im the figure. Input data is included in Appendix A.14. 


Pore pressures obtained by ADINAT using the above cited 
theory are shown in Fig. 5.11. Finite strain analysis is 
also performed and compared in the figure. Again, there is 
little difference between these results. This same problem 
is also solved by Raymond (1969) with small strain 
approximation and good coresspondence in the solution is 


noted. 
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material laws 
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Figure 5.11 Pore pressure distribution in thin layer analysis 
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5.6 Poskitt Analysis 


Objective: 


To illustrate ADINAT application to Poskitt's (1969) 


analysis. 


ADINAT Nonlinear Equation: B10 


Analysis and Results: 


Poskitt has given the same formulation as equation B10 
in Fig. 4.4. In order to make use of the given material 
parameters in his paper, changes of variables have been made 


to equation B10 as follows. Setting 


c = (ma)/h (5.9) 


and 
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and substituting these into equation B10, yields: 


m 
or Carp) Ou Ou 
kh eee ke ea Py (Smit) 
CCue a (i -Les) Oc OT 
f 


ADINAT analysis is carried out with the parameters a, 
precee CVi7n- = Ofs?2, 0.37, 0205, 15535x10- */mim respectivelys 
f 
Input data is in Appendix A.15. The output pore pressure at 


T=Oe5 1S shown in Fig. 5.12. 
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Figure 5.12 Pore pressure isochrones in finite strain - thin 


layer analysis 
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5.7 Gibson-England-Hussey General Formulation 


Ob ject.ive: 


TO verity the application of ADINAT with Gibson ef ai 


(1967) formulation. 


ADINAT Nonlinear Equation: A7 


Brief: 


These authors have given a general treatment of one 
dimensional consolidation of a single layer soil. They 
arrive at a governing equation for a thin layer which is 
essentially identical to equation A7, but the variable 


ecoefiacizentsc ve 
F 


ny Conc) 


is assumed either aS a constant, or aS a linear Function oF 


the, void ratio e, 


| Sie 
‘agcy Lexvensp » Meuip aved avodguepeeen 


pm 
be 
‘“ 


a4 ;m i 
a P i 
ae na ‘tavel olonia:s/28 @mektaedi loanc 

xf nic 6 to1 Aoltaepe Pal 
| : 5 ae 
— _ a _ a — ‘ r -— 
7 a" ‘h. nobigapeet.iaotinebi qital 


c SIG G7 + m(e-eg) (5213) 
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Co and m are constants which constitute the mode of 


VarlativonVo£CC 
F 


Description and Analysis 


For comparing the results of Gibson et a] and to use 


their given parameters, equation A7 is tranformed into 
0 de! de' 
—|c: — | -— (5.14) 


by changes of variables, where a' and e' are ratios of the 


orginal values. T' and C' are defined below, 
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Thescase of A=0.4 and e /eog=0.2 with the linear 
se 


material model, Eq. 5.13, is illustrated here. The datafile 


is given in Appendix A.16. 


Results: 


ADINAT predictions for the particular case taken from 
Gibson et a] is presented in Fig. 5.13. Only half of the 
layer is shown because of symmetry. Gibson et a/'s results 
are also reproduced in the figure and good correspondence is 


found. 
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5.8 Discussion 


It can be seen from the preceding examples that most 
thin layer theories can be retrieved from the ADINAT 


equations in Fig. 4.3 and Fig. 4.4. 


For cases where analytical solutions are obtained, 
excellent agreement can be observed with ADINAT predictions. 
For others in which the original work has been further 
approximated in some way, good correspondence is also shown 


between the theories and ADINAT results. 


The closeness of the results of the finite strain and 
small strain analyses may be due to the insignificance of 
strain or void ratio changes in a thin layer problem. An 
exception to this occurs in the example of Davis and Raymond 
(1965), where the greater discrepancy is probably due to the 


larger load ratio used. 


The difference between the finite strain and small 
Strain analyses becomes more notable in deeper layers. This 


will be clearly seen in the next chapter. 
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6. NONLINEAR CONSOLIDATION ANALYSIS OF THICK SOIL LAYER BY 


ADINAT 


6.1 General 


In the following, examples are given for cases of thick 
layer analyses. Only a few papers on this subject can be 
found, probably because of the necessity of dealing with 
highly nonlinear equations. General descriptions, including 
the finite element model, given in Section 5.1, are also 


relevant here. 


Notice that eyo iS not a constant in thick layer 
analysis and the initial effective stress in the Lagrangian 


coordinate requires 
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The effective stress and excess pore pressure relationship 


at time t are given by Eq. 4.35 or Eq. 4.42, depending on 


the coordinate system used. 


By integrating Eq. 4.40 and assuming the logarithmic 
effective stress - void ratio relationship, the following is 


obtained: 


o> Ce 

a = (1+teg)b - ——— (1 - ¢ + clnc) (653, 
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W 
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eae 

Gea + *b (6.4) 
O¢ 


The above expression shows the relation between the 
coordinate a and b, and it is useful in the thick layer 
analysis with the logarithmic compressibility law. Notice 
that the Lagrangian coordinate @ in Eq. 6.3 can be 
substitued by the Eulerian coordinate Z, in which case the 
Same change also applies to Eq. 4.40 for computing the 


reduce coordinate D. 


Several other symbols are employed in thick layer 


analysis. The initial void ratio at surface, eg, can be 
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The surface load ratio is computed from 


and the self-weight ratio § is defined as 
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where £ can be obtained from 


© represents either the coordinate @ or Z. 


using the average submerged weight, termed the average 


self-weight ratio, 1S expressed as 
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SR=oeChyeeht)lo/ oc; (6.9) 


in which y is the average submerged unit weight of the 
b 
soil. 


In applying thick layer formulations to ADINAT, the 
compressible stratum is divided into contiguous layers of 
finite depth. The expressions for material parameters are 
then evaluated at the mid-point of each layer. Although 
this is only an approximate method of the application of the 
thick layer formulation, good results are obtained as is 


Shown in the following examples. 
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6.2 Small Strain Approximation 


6.2.1 Raymond Analysis 


Objective: 


To illustrate the use of ADINAT in Raymond's (1969) 


analysis. 


ADINAT Nonlinear Equation: A4 


Brief: 


Raymond has studied the consolidation of deep deposits 
in the context of small strain. Variations of permeability 
and compressibility are taken into account by employing 
logarithmic relationships. Notice that Raymond's 
dimensionless equation Eq. 11 used in his calculations is a 


further approximation of his general equation Eq. 6. 


The formulation Raymond has given to solve deep deposit 
problems may be considered as a pseudo-thick layer analysis. 
The variable term (1+te,) is incorporated in the time factor 
to whiche theseoutputetime!l for poret pressures results are also 
referred. This is of little practical use because the real 


time may not be easily determined. 
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Description and Analysis: 


An example is shown with the same variables assumed by 


Raymond. Setting 


Ko, lm lO(1te,) 


T = —————_ ¢t CGraudy 
yooh? 
wW 


and with the logarithmic material laws, gives from equation 


A4, 


— [ oa! — J = —— — (6.11) 
One pes a! oT 


where z' is the depth ratio z/h. The soil profile and 
material parameters for ADINAT analysis are given in 


Evgssoele. The datatile is contained in Appendix JA. i7. 


Results: 


Pore pressure distribution at 40%, 60% and 80% 
consolidation computed by ADINAT are compared with Raymond's 
Solutions in Fig. 6.2. Slightly larger differences can be 


seen at later consolidation stages. 
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Figure 6.1 Nonlinear consolidation analysis in Raymond (1969) 


example 
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Figure 6.2 Pore pressure distribution in small strain - thick 


layer analysis 
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6.2.2 Davis Theory 


Objective: 


Toddemonstrate the lise of ADINAT imeDavius (1971) 


analysis. 


ADINAT Nonlinear Equation: A4 


Brief: 


Davis has given a small strain theory for a thick 
layer. His governing equation is essentially the same as 
equation AG except that the coefficient of consolidation is 
assumed constant and the coefficient of volume 
compressibility is inversely proportional to the effective 
Stress. Closed-form solutions are provided in some limiting 


cases along with these simplifications. 


Analysis and Results: 


Substituting the above mentioned approximations to 


equation A4, which becomes 
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— UG 12) 


This equation is then used for the analysis. An example of 
the load ratio, a, of 1.005 and the average self-weight 
ratio, §, of 5 is given in Fig. 6.3. The input datafile is 
in Appendix A.18. Pore pressure results at time factors 0.2 
and 0.3 given by ADINAT are plotted in Fig. 6.4. Davis' 
closed-form solution for a, = 1 and Terzaghi's solutions are 


also compared in the figure. 
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FINITE ELEMENT MODEL 
FOR THE CLAY LAYER 


MATERIAL MODEL: 


: NODE 
e - log a 
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CONSTANT Cv 

zZ 

MATERIAL PROPERTIES: 

Cv = 1.0 m?/day 

yY = 2.0 t/m 

b 21 


AREA = 1 m? 


Figure 6.3 Nonlinear consolidation analysis with Davis (1971) 


formulation 
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Figure 6.4 Pore pressure variation in thick layer analysis 
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6.2.3 Viggiani Numerical Solution 


Objective: 


To compare ADINAT predictions with Viggiani (1973) 


numerical solutions. 


ADINAT Nonlinear Equation: A4 


Brief: 


This author has given a numerical treatment to Davis' 
(1971) formulation, and also considered the effects of 


variable load increment with depth. 


Analysis and Results: 


The soil layer and material parameters are shown in 
Pig 0.2.8 be darainle is in Appendix Al19.0%8q. 6.12 is 
used for ADINAT analysis and pore pressure results at time 
factor 0 27and 0. 3yare given in Fig. 6.6.) Wiggrani's» and 
Terzaghi's solutions are included in the figure for 


comparison. 
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FINITE ELEMENT MODEL 
FOR THE CLAY LAYER 


MATERIAL MODEL: 
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Figure 6.5 Example of nonlinear consolidation with Viggiani 


(1973) data 
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Figure 6.6 Pore pressure variation in thick layer analysis 
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6.3 General Analysis 


6.3.1 Example with Logarithmic Material Law and Constant 


Coefficient of Consolidation 


Objective: 


To illustrate ADINAT application to nonlinear 
consolidation with logarithmic effective stress - void ratio 


relationship and constant coefficient of consolidation. 


ADINAT Nonlinear Equation: A4, A8, AQ 


Analysis and Results: 


The problem layout and material parameters are 
presented in Fig. 6.7. The small strain equation Ad)and 
both of the: two finate strain - thick layer formulations, AS 
and AQ, each of which utilizes a different coordinate 
System, are used for analyses. Input data can be found in 


Appendix A.20. 


Pore pressure distributions predicted by ADINAT are 
Shown in Fig. 6.8. Greater difference between the small 
Strain and finite strain formulation can now be seen. 


Terzaghi's solution is also included in the figure. 
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MATERIAL MODEL: FINITE ELEMENT MODEL 


FOR THE CLAY LAYER 
eq- log ca 


CONSTANT Cv NODE 


Tole m 


MATERIAL PARAMETERS: 


a, b 
@o = 0.9% Cc’ = 0.4 
Cv = 1.0 m?/day 
Yao 42.0 0/7 mF 
Yoo ge 1.0 /m 
Y¥, = 2.9 gt/n 21 


AREA = 1 m? 
REDUCED CLAY LAYER THICKNESS &¢ = 8.438 m 


Figure 6.7 Consolidation analysis with logarithmic material 


law and constant coefficient of consolidation 
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Figure 6.8 Pore pressure distribution in thick layer 


analysis 
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6.3.2 Gibson-Schiffman-Cargill General Analysis 


Objective: 


To illustrate the use of ADINAT in finite strain - 
thick layer analysis with Gibson et a/J's (1981) general 


theory. 


ADINAT Nonlinear Equation: A4, A8, Ag 


Brief: 


These authors applied the same general equation, 
Eq. 13, from an earlier paper, Gibson et a] (1967), but the 
eneectsoieself-weaghterst includedeanvtheirai98 Ispapers 


Their governing equation 


oe Me Sere ,) k do’ de de 

——_ — —]+—= 0 (6813) 
y de 1+e db db y (ite) de ab dt 
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(Site 


is basically the same as equation AQ. They made further 


assumptions and arrive at an approximate equation as 
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my. y ) == (6.14) 


where A' and g are assumed constants, which may be computed 


as: 


Bae 1n[——————_] (Gren: 


quamey / [i+tela)] Coral 


The bracket (a) denotes the average value and (~) refers to 


the final value at the layer bottom. 


Description and Analysis: 


The soil layout and material parameters from Gibson 
et al's (1981) example are given in Fig. 6.9. The 


compressibility law employed here is the same as that in the 


paper: 
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Figure 6.9 One dimensional nonlinear consolidation with Gibson 
et al (1981) analysis 
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The coefficient of volume compressibility can be determined 


Wiuitethis melation and /Eq.\3 J8, 


m = A'[Le-e(~)]/(i1t+e) (6.18) 
Vv 


Note that Cv is treated aS a constant in the paper and hence 
the coefficient of permeability k can be calculated from 


ie &S OK Pa ae 


Twenty solution time steps are used in this example. 


Datafiles are given in Appendix A.21. 


Results: 


Pore pressures predicted by ADINAT at two different 
times are shown in Fig. 6.10 and Fig. 6.11 with Gibson 
et al's (1981) solutions. Reasonable agreement can be seen. 
A small strain analysis of the same problem is given in the 
figure. Terzaghi's solutions is also included for 


comparisons. 
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Figure 6.10 Pore pressure isochrones in thick layer analysis 
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Figure 6.11 Pore pressure isochrones in thick layer analysis 
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It can be seen from the figures that the conventional 
theory is an underestimate while the small strain 
approximation is an overestimate of the pore pressure. 
Little 'bumps' on the top and bottom parts of the smooth 


curves are due to the limitations of the plotting machine. 
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6.4 Discussion 


Since thick layer analysis allows for variation of 
effective stress with depth, and large changes of void ratio 
Or strain, the difference between small strain and finite 
Strain analysis becomes more significant as shown in the 
preceding examples. In addition, pore pressure isoschrones 
are all skewed because of the unsymmetrical stress 


distribution in the soil layer. 


Note that the average submerged weight is used in the 
alalyses by Raymond (1969), Davis (1971) and Viggiani 
(1973). Since this value does not consider the variation of 
void ratio throughtout the soil depth, it is apparently 
inadequate to take account of the self-weight effects. The 
theories derived by the cited authors, therefore, are only 


approximate to the thick layer analyses. 


In the analyses of thick deposits using ADINAT, 
preparation of input data may become tedious since material 
parameters change with depth. Simple supplemental programs 
can be eaSily developed for this purpose. Examples are 


given in Appendix B. 
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7. CONCLUSIONS 


Application of ADINAT to various one dimensional 
consolidation theories are illustrated. Within the scope of 
the present work, finite element solutions are readily found 
with simple procedures for linear and nonlinear analyses. 
The overall simplicity in the use of the program makes it a 
convenient tool for general consolidation studies. Although 
examples are selected from a large area of ADINAT 
applicability, only part of the many possible applications 


of the program have been shown. 


The high accuracy of ADINAT and its incremental 
solution procedure are particularly useful in soil problems 
with nonlinear material behaviour. The ADINAT equations in 
Fig. 4.3 and Fig. 4.4 provide a simple means for extensive 


research into nonlinear consolidation behaviour. 


Use of the program for nonlinear analysis is 
Straightforward once the material relationships for a given 
problem is determined. These relationships are then input 
into the program at discrete points in the range of 
consolidation. Linear interpolation is carried out by 


ADINAT to obtain values between the input points. 
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If the effective stress is a function of void ratio 
only, there should be no great difficulty in finding their 
functional relationship experimentally. A series of 
Standard consolidation tests may well be performed to 
accomplish this. The problem of determining permeability 
and void ratio relationship appears to be a little more 
difficult. This may possibly be resolved by additional 


measurements attached to the consolidation equipment. 


Parametric studies can be undertaken easily by the 
Program. Eifects of load ratio, self-weight ratio, soil 
depth and other physical conditions can be investigated. 
The assistance of a small supplemental program for the 
preparation of input may greatly simplify the use of ADINAT 


for this purpose. 


New sets of useful consolidation parameters for 
nonlinear theories may also be developed by reviewing the 
proposed formulations. Comprehensive experimental and case 
Studies, however, should be performed to demonstrate the 


usefulness of these parameters. 


The operation of ADINAT can be stopped and restarted at 
preselected time steps to allow for changes of input 
options. This restart capability may be used for analysis 


of problems with stress history and time dependent loading. 
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The program also permits the elements to be active or 
inactive at any specific solution time. This option may be 
used to simulate the physical situations for a construction 


period or moving boundaries. 


The orthotropic material model of ADINAT can be used to 
analyse anisotropic conditions in soils. Variable load 
increment can be handled easily by the program by specifying 
the appropriate loading at nodal points throughout the 


layer. 


Impeded drainage boundary conditions may be treated as 
layered problems with the additional top and bottom layers 
being assigned appropriate permeability values. Radial 
drainage conditions can be modeled using axisymmetric 


elements. 


Assessment of the numerical behaviour of ADINAT, such 
as stability in nonlinear analyses, is required in future 
work. A better understanding of the numerical behaviour is 
essential for further development of the program. This will 


also facilitate program use and reduce operating costs. 
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The ADINAT program (a Finite Element Program for 
Automatic Dynamic Incremental Nonlinear Analysis of 
Temperatures) is available on the MTS system at the 
University of Alberta. A copy of the user's manual can be 
obtained from the Department of Civil Engineering. The 


following commands can be used to run the program 


$RUN 3012:ADINAT 5=input 6=output 


Listings of datafiles corresponding to the examples 
discussed in the thesis can be found in this appendix. The 
first card in each datafile contains a brief note of the 
file. Computer programs used for the preparation of input 


data of several examples are given in Appendix B. 
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A.8 Linear variation of modulus of elasticity 
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A.9 Sinusoidal variation of compressibility 


1 


TZBDASETESOIZIASETBSOIZIASETEPOIZIASETESOIZIASETESOIZIASETEIOIZIASETBE9O1Z34S67880 


1-0 CONSOLIDATION 


4) 1 
1 ° 
10 

4316 
1 1 
4 ° 
ao ° 
ai 1 
1 1 
1 

1 

1 

1 

1 

1 

1 
1 2 
1 2 
° 

1 4o 


964666166 
2 
ns 
6984216342 
3 
s) 
824418710 
4 
1 
785703895 
S 
1 
-€ 88454649 
6 
a] 
623206236 
7 
1 
$602408283 
é 
4 
.@49998E6 790 
s 
1 
442815446 
10 
1 
389078329 
11 
1 
.3381098742 
12 
1 
293214762 
ves 
a] 
.251677349 
14 
1 
214753596 
1s 
1 
182671146 
16 
1 
185627798 
5 es 
1 
133780287 
18 
a 
117293248 
18 
1 
YOEZIEIES 
20 
if 
-1O0B938E7 
au 
1 
100693867 
a2 
a] 
- 106232389 
23 
1 
1197293248 
24 
1 


te) 


EO 
EO 


=e O 
Be ed 


EO 


A Ke) 


° 


ooo°o 


cok uM NW 


2 3 


1 1 


oooo 


EO 
5A ae) 
EO 
EO 
.EC 
EO 
.EO 


Le) 


(SIN VAR My, 


FP nee 


° 


ob DOR) 


313) 
.EO 
. EO 
at Shey 


EO 
eC 
Ec 


4 


10 


19.5 


wUM BER 


° 


Gb) ae es an eed 


EO 


.EO 


EO 
EO 
EO 


BEC 


5 


ee 


6 


EO 
EO 
EO 
EO 
EO 
EO 


ao 


a 


oT 


eal he) 


EO 
EO 
EO 


.EO 


Ec 


DATA:Schif ¢manéGibsonéss 


208 


“waif tatzeas qmo 


Te 
oe ; 4 s . as 
* a 
{ Sea anes rbtataungeyseecseneesnants ere 
es ayaa Twn erees eee gee a ad few dias 
; + soguaenenamntts ci eien™ ¢ es 
: : 


- 


rie 


23.t *? | ad - et as.¢ uv. 
, es? sa oé.¢ i a eI 
. oa f abt ofr. i of 
39 os oe. ei. 45 1 os 
ee 5 es .f eg 7: “2.4 o> 
» os 6 as? oe a on. 23 
esr sa ee 
: = 8 
r) . 2 s 
©. 4 . .- 
? 
- ' 
¥ Q r , 
cnasinvake 7 
’ 4 7 “7 : 
' a 
erthe ance 
, , & 
A 
peorseset 7 
- i - e dss 
> : i | 
salt. feuscenaa 
: * 
| Pa veeuegets : 
i. % z ae 
¥ 7 
tieseneda. 4 
éovaséeus vl 
.:.,.. = 
ad 
easyacan aa 
‘ = 
eS i 
be iecuvedee oat 
| ' ? 
? 
Ss tees ccs 
, = 
cersicea 
; : . 
i . 
etriers€ oo) 
Fr *: 
in 
deaeetere: 
: ’ ar 
4 oe 
savrevater 
’ a 
; seretagas 
ae Oe 
rd vengercer 
d ; - ’ _ . a 


pa 


CARD 
NUMBER 


.133790287 
25 

az 

.155627798 
26 

ihe 

.182671144 
27 

ihe 

.214753594 
28 

Ab 

.251677343 
29 

te 

.293214762 
30 

the 

.339109742 
31 

1 

.389079329 
32 

Be 

442815446 
33 

Ny. 

.493986790 
34 

he 

.§60240883 
35 

fi: 

.623206236 
36 

; ae 

.6884946439 
37 

tis: 

.755703595 
38 

ee 

.824418710 
39 

Alp 

.894216342 
40 

A; 

964666166 


word nanaWwn — 


10 


12345678901234567890123456789012345678901234567890123456789012345678901234567890 


OnMITDRW PUN — 


10 


OnIMH PWN 


10 


WOeOIMMH BWN = 


10 


c.0 t UM N 


No UUM! B87 ER 


209 


Erna 


* 
- 


ror 
: : Dies 
on ' 


erp se abinagiassipcr=*~s said sg rAd a pl wee Bi 


ee — S} case rang ease 


> a 


qeessk: ba. - bee Eaenn §.. 
ee a. 


: SELTaaeMRAAES 


I 
- 


32 


- 


‘7 


Ws @@'s 


b = 


=e ee 


‘ Wes cei 
areas « 


steapabe? iat 


' 
+ (deneeds ae 
«= % 


A. 


7 
edeeentn ; 
“ 

' ve 4 4 
ore : 
iW oe i 

‘s ; 

i 

—— 

« 

‘oF 
‘sors ey 
a 

a ecu 

#7 

6 

a ve ;* 


A.10 Constant coefficient of consolidation with variable 
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A.15 Finite strain-thin layer: Poskitt (1969) 


cok UM W WUM BER 
1 2 3 4 5 6 / 8 
TZI4SETESOIZIASETESOIZIASETESOIZIS4SETESOI1ZIASETE9OI23ZASETE9OIZI4SETEHOC1ZI“ZSETESO 


1-D CONSOLON (FIWITE STR, THIN LYR, D.DR) DATA Poskitt6S EO: B10 


21 ° 1 1 a] ° 10 
1 ° 
10 
.Os8 
1 ° ° ©. ° 1 
21 ° © © 3.1615$265 
° 1 
1 2 ° 2 
1 2 
° 1 Le) 1 
1 1 1 ° ° 
21 1 1 ° 
1 20 1 2 1 42 ° 
1 1 
oo © os o 10 ©o.18 © 20 oOo 25 ©. 30 o. 35 
oOo 40 o 45 © $0 © $5 © 60 o 65 o. 70 OE7s 
© 80 o 8s © $0 °o.95 1 002 467401102 436756182 40651083 


") 


3766541462 JB471E2ETZ 318O091S64F BZUSIIMESE? 2610297463 2330494572 205425772 I1781ES9E6 


2 181282322 1246687262 OF 8463182 072576642 047023172 02179839) 89688798) $7231770 
1.948053331 92410073 
1 
1 1 2 1 1 ° 
20 20 2) 1 ° ° 
STOP 
1 2 3 4 cd) 6 7 8 


TZIASETESOIZISSETESCIZIASETBIVOIZIGSETESOIZIASETESOIZZIASETESOCIZIGSETESOIZIESETESO 
CO ak 1) AN NEU 4 78 ER 


Leo 


ra 


i, 


ee. a ne . ae i 
SV tus cr eer 
ease yeyentbentarans) begrenseet eae 


ee See te Dnt beta haliteliehae ee 


3 ea 4 


ane te eee, Te ‘a tot 7 o ya 0 sia 
eto eno on @) Of OP 8 ae © oth E aa > *~*s - “Seo, 
cye( seem ra) eervete Bet) Osrss t20.8 pe 2enme za @ o a a 
azovracy ctheseeun S2eesor’s Dat eeei Se THESE s6tOtoy ots car Pow Ee 


raeeareeo srrosetee ba as ” Bkse8a 5 ecngt €eCcerr 


ee ae <2 1ctoree er. 
T4o1e 


iS - "ysis 


ey 


ee eee ew aaa ae vinse 


ee ee ee ee 


. * a ¥ s 7 -¢ s * r 
of iepenic manner Nana Wiser «NR eae EE 
a ee ee ~ ? 2a } A 

iy a 


226 


A.16. Finite strain-thin layer: Gibson et al..(1967) 
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Small strain-thick layer: Raymond (1969) 
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1-D CONSOLDN (Yerzaghi’s Analysis, D.DR) DATA:Davis71(Tv=.2) 
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1-D CONSOLON (Terzaghi’s Analysis, D.DR) DATA: Viggiani73(Tv=.2) 
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1 1 ° loan °°. ° 
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10. 10. 10. 10. 10. 10. 
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A.21 Thick layer: Gibson et al (1981) 
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099035482. 
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128658202. 


3.00 
Ea Oe 


<2eaeo256=5 ; 
Q9168O06GE-S. 
2ZS4960E-4. 


3.00 
7.00 


-008350853. 


031511852 


118909624. 


3.00 
7.00 


-237880E-S. 
8S2074E-5. 
2E68460E-4. 


3.00 
7.00 


007718100. 
029124163. 
109899706. 


3.00 
7.00 


Z21964S9E-S5. 
7T91S88E-S. 
-252590E-4. 


3.00 
7,00 
007133291 


026917392. 


101572480 


3.50 
7.80 
3799 73E~5 
132243E-4 
377613E-4 


= feat Ko) 
7.80 
012487778 
047122491 
177816193 


3.50 
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3S51981E-5 
123228E-4 
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3,50 
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164342858 
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7.80 
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3.50 
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3.50 
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321643E-4 


3.50 
7. 50 
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140381503 


3.50 
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009111781 
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12989744637 


3.50 
7.50 
258748E-5 
S239 75E-5 
Z867S7E-4 


3.50 
NiO 
-008421371 
031777949 
-119913738 
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1-D CONSOLON (Terzaghi’s Analysis, D.DR) DATA:Gibson-etal8s!1 


aut 1 ° 1 1 }. 20 
1 ° 
20 
180.9 
1 1 °. ° ° ° 
2 ° °. ° 2 1 
21 1 °. ° 10. 
1 1 
10. DOs 10%, 10. 10. 10. 
10. 10. Lg oaks 10. 10. 10. 
TOT, 10 10. 10. 10. On 
10. 
1 2 ° ° 
1 Fe 
©. °. SOOO. °. 
1 20 ° 1 1 ° ° 
1 1 
-001944 
is 
1 1 1 1 ° 
20 20 21 1 ° ° 
STOP 
1 2 3 4 5 6 Ci 8 
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Computer Programs 
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LIST OF PROGRAMS 


Program Page 
B. | pemgechis classical analysises. te ceed rere aera oh!) 
Ba 2 CEeaymogtuWwOn layered analy Sic. x nice see te sie etre mere 260 
Bees Maoh SimiOLOSed-fonmesolution a... .s<:6 sss ee1e ele eee 261 
B.4 HGcuMaN iw eed OSEA-fOrm "GOlUEs OM weter mee cre eer reie ee 263 
B.S Dota ne le he U2 EC. GAGs i a sre eile esr eens ere eee te 266 
B.6 Datel Met Ae OS" EC 6 AG etoeres theme ener at eek name ee ewes eee 260, 
Bas Lena Lind Cerra 20s IsCh eos ptereBele ot aiete saclerenwretetoectaom Cad ereners 268 
B.8 Da ba feu ei Ags lets cECh ow AG sea tepevewtdamn nhs recBinc 34.0 elo haber elenatelene 269 
Bao Data fislem Acc lis, EC. ACs scslcheteee tel imino teaei cece ee i Pose generac 270 
Baa Dertaret Gaia c ks Es. WAS reravel ee tenetereconet sn menebarene cr cures Os 


Note: Most of these computer programs are 
written for the illustrated examples in this text. 
Modifications of these programs are necessary for 


their general use. 
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B.1 Terzaghi's classical analysis 


SSSSCSSCSSSCSTSTSSTSRRESTSSSSESESSSSSESseSeeseseseuesssesessaees 


s 
see TERZIAGHI2S see = 
This program is to compute the excess pore pressure s 

et various depth and time aeccoraing to Terzaghi's s 
s 

s 

s 


classical theory of consolidation 


SSesceaeeseesesesesreesesseesessessesseeeessessesessessesessseaeesese 
COUBLE PRECISION 2H(28),TV¥(16),UE(15,26), aM, AMT, AMZ 
DOUBLE PRECISION DSIN,DEXP,DFLOAT 

Tv (t)80.0800 

0O 8 Ive2,4 

TV( TJ) ETV(1U-1)4%0. 0800 

00 10 I1J38$,11 

Tv¥(T1T)@eTV(11-1)46.100 

2H(1)20.08D0 

0O 168 JJ#2,20 

ZH( JJ) EZH(II-1)90. O08 DO 

OO 20 181,11 

00 20 J#1,20 

VE(I,J)€0.00 

BDO 20 Ke1,60 

KMeK-1 

AM= (2. DOSDFLOAT(KM)+1.00)#3.141692688400/2 DO 
AMZEAMEZH( JU) 

AMTs-1 .DOFTV(])shamMa2? 
VE(1,J)SUE(I1,J5)*DSI1N(AM2)*DEXKP(AMT)82.0D0/am 
CONTINUE 

WRITE(6,60) 

WRITE(E,61)(TV(I),121,11) 

OO 25 L2#1,20 

WRITE(G,62)Z7H(L), (UE(T,L),1#1,11) 
FORMAT(’'O’/’ Deptn’ ,464xX,’'Time Factor Tv’) 


FORMAT(‘ Ratio’ ,11F9.2//) 
FORMAT(FPE.2,2K,11F9 6) 
sToeP 

Exo 


B.2 Gray's two-layered analysis 


BPRASTESSLCKNATASseseseeeseseeseeseseseseseseseeeeosseeeesseesteseseeseseesaesesseseses 


e s 
* ees Grays& #82 s 
. This program is to find the analytical solution of s 
s tne excess pore weter pressure for two-leyered soils s 
s according to Gray H. (1846) for certain values of sigma, s 
s mu, lembde as given below s 
s z 
s a 


DOUBLE PRECISION AW/(10),TV(B),2m4(28) ,VE(S,28),CWH,PARC, PARMO, TPP 
DOUBLE PRECISION PARAZ,AMU,S1GMA, ALAMDA 
AW(1)80.0714326347 

AW(2)8®O 215 

Amw(3)80 3604 

Aw(4)80 5078 

aAw(S)2#0.6571 

Aw(6)80 8076 

Aw(7)#0.96533 

Amw(8)2£1.1157 

S1GmMAt2 © 

ALAMDSr4 © 

AMUz5 © 

Tv(1)#20 0 

Tv(2)#2100 0 

T¥(3)#200.0 

ZH(1)86 06 

Bo 160 1182,21 

ZH( TT )e2M(1y-+1)°0.28 

po 20 181,3 

BDO 20 v#1,21 

UE(1,J)#0.0 

oO 20 Ke1,8 

PARCEAMUSALAMDAS AN (K) 

CNE2 SDCOS(AN(K))/(SIGMASAN(K)SOSIM( PARC )822¢PARC#DCOS(AKiM) ) #82) 
THPPE-1 8T¥V(1) SAN (K) SANK) 

3F- GS. GE..18)) 69, TO 16 

PARMATAMUSAN(K)8( 1 *A2LAMDA-TH( I)? 

VEC I, J) EVE( IT, J)eCHSDEXP( TNPP HF DCHS(BwWIK ) LSDSI MO ARME) 
GO TO 20 

PARATSAWN(K)8ZK( J) 

VE( 1], J) SUE( 1, J)eCM2DEXKP( TNPP) SP DCOS(PARAZ)SDSIN( PARC) 
COWNTIMUE 

WRITE(E, 83) 

WRITE(6, 60) 

WRITE(E G61)(TV(1),T8#1,3) 

Do 28. tat 21 

WRITE(6,.62)Z2HM(L),(UE(T,L),181,3) 


sTOP 
FORMAT(’ Deptn’,8X,’Time Factor Fw! 
FORMAT(’ Ratio’,3F8 O//) 


FORMAT(F5 2,2xK,3F9 6) 
FORMAT(‘'O'//4K, Analytica) Solution by Grey’//) 
END 
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B.3 Martins" closed-form solution 


SSTESeSSSeSErseseesesetc sess sereesrseesseeeseseersetsreasegeseesesazaeseeegsssss 
zat MARTINSES #28 


* One dimensional consolication with linear spatial 
variation of permeability # 


This program is to calculate the excess pore pressure 
at various cepth eat time t according to Martins’ 
analytical solution im 1965 The soil leyer is assumed 
to drain freely at both top and pottom. 


s 

z 

s 

= 

z 

= 

s 

s 

z 

2 

= 
Input variables 2 
T = Time required, sec 2 
KO = Tre coefficient of permeability, cm/sec bd 
MV = The coefficient of volume compressibility, cmt#2/g 2 
BETA = A dimensionless contant relates to the variation = 
of permeadility by the eqn: K © KO(1+#BETA#7Z) . 

GAMMA = Unit werght of soit, g/comst?3 » 
GAMMOW = Unit weight of water, g/cmt2s3 z 
H2 = Vota) Chirckness ‘of Sou), cm s 
vo = Initial uniform pore pressure, g/cmtt2 2 
NOEL = WNumber of equa) segments of the soil profile z 
specifrieo for the outputs s 

AN(n) = Roots of the eqn. JO(A) YO(RKA)-JYO(KA) YO(A)=0 s 
where k is the variable LAMBDA in the program s 

NORT = Total number of the above roots used s 
z 

= 

7 

= 

z 

z 

s 

s 

z 

s 

s 

s 

2 


Use of the program 

1 First to compilde the program Dy: 

RUN *®FORTGTEST SCARDS=MARTINSES SPRINT=-1 PAR=10,SOURCE 
2 Repeat the following for subsequent executions 

OEBUGC -LOADH+s IMSLDPLIB S=datafile 6=output 


Remarks: 

Tne routines MMBSJO and MMBSYN in the IMSL Library will 
be wsed Tneir descriptions are shown at the back of this 
program 


re ee ee 


BUFAKRTSERTFTSRKKC RESTS eeeeerseeserrssaeseeseeaeaerseaeesueessesaseseses 


IMPLICIT REAL*8(A-K,O-7) 
REALL*®&6 KO,MV,LAMBDA ,OSORT,DEXP ,DFLOAT ,MMBSJO 
DIMENSION AN( 10) ,ALA(10);,BY0LA(10),UT(10),21100),2H( 100) ,UE( 100) 


zseeeepeazazreseseseeeseesetreseesaeseesvseeseeezsesreaeeaeeeeeeeseeeeaesetrxreesrseerese 
z z 
* Read 81)1 input deta z 
z = 
Basaaeeezsrsauezttezreresetaseereeerezpeesereraresesveererexseseaueeeeecezresereee 
READ(S SOT, KO,MV,BETA,H2 
READ(S S1)*CEL,NORT,UOC,CAMMSL ,_GCAMMAW 
READ(S 52) (4N41),}=13,NORT) 
Seseererereeseseexeeeteseeeescesateaeresesreeeveeeresrsreseestreereeseeeeert 
s z 
= Calculation of the average consolidation at the given 2 
2 time T s 
z s 
Bezpezrersrsxarsreeeearseeseeeeseeetreesaesetreresetreeceaeaeseeasereerseeerees 
LAMBDA=DSORT(1,.DO+BETL) 
CHI=BETA®BETA#KO/(4 DO®H2*=H2* GCAMMAW?FMY ) 
ACT=0 DO 
po 15 J=1,NORT 
ALA’‘1T+ = LAMBDA ANI I) 
BJOA=MMBSJO(AN( I), ,IJEJ1) 
BYUOLA(} )EMMBSUCLALA!T*, 1EJ2) 
EV=DEXP(-1 DOFCHIFAN(I)®AN(1/#T)/(BJOL+BIYOLA(I)) 
UT( 1] )SBJOAtTEYU 
ACTZ=ACT+(BUOA-BYOLA!)) 12 EU/t(AN( 112 AN'I1)) 

15 CONTINUE 
LYVCON=4 DOFACT/BETS 
seeressxrsseraeeraereerereneraeeererereereaerreerereesrereresrereaeeaeeeE 
5 . 
2 Compute the excess pore pressure at various deptn a 
= s 
Sxesseeesrsrrsesssresezrseseesseesesr sr xezareseereseexeEwsrsrsreesresssrsserszrisese 


Pl=3 TOISSZ2ESBSSSTSIZISAEZEAIZBEZTZTIGSOZSBSHISTIEIS3ZSS37TS5IING 
NONOD=NOEL+1 

DIV=H2/DFLOATINOEL) 

2 1 v=. DO 


00 20 J=1,NONOD 

Zi dt IEZ(JI+DIV 
ZH(J)=EZ145)/H2 

R=DSORT(1 DO*BETAFZHI VU) ) 
UET=0O DO 


po 18 L=1,NORT 

RA=R*ON'L) 

BJORD=MMBSVO(RA,IEJI3) 

CALL MMESYN(RG,O DO,:,YORA,IEY!? 

CALL MMBSYN(ALA‘L:,©.00,1,YOL4,IEY2!) 

UETZUETS#UT(L)® (BYUORDL* VYOLA-BYOLA(L)* YORA) 
18 CONTINUE 
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VE(J)=U0O*#PI*UET 
CONTINUE 


RRR EKA KK SKK REA EAR ERAT EAE EERE EERE HERE EREEKEREEEKEKEET 
> = 
* Output quantities * 
x 

+ 


HERE RAK ESET H ERE AERA A TAR HAKEEE EKER EEREEEERERKEEEKES 


WRITE(6,60) 

WRITE(6,61)AVCON 

WRITE(6,62)KO,BETA,MV,H2,GAMMA, GAMMAW, UO, LAMBDA, CHI, NORT, NOEL 
WRITE(6,63)T,(Z2(1),Z2H(1),1,UE(1I),1=1,NONOD) 
WRITE(6,64)1TEJ1,1EUJ2, 1EU3, IEY1, IEY2 


STOP 


SRST E REESE AERA EERE EAEE EEA KA AREA ETH EE HEE TAR EKEKETETAEXET 
= 3 
* Format statments * 
= x 
RRR TERE K KSEE KAKA HEE ESHA KSEE A KEK AK KAKKKK KEE ERE EE EEKREKREEEEXKE 


FORMAT(5G20.0) 
FORMAT(215,3G620.0) 
FORMAT (3F20.0) 


FORMAT(‘1°////6X’The One Dimensional Consolidation with ° 
i’Linear Spatial’/18X’Variation of Permeability’ 
1//13K’Analytical Solution by Martins (1965)‘///) 


’ 


FORMAT(4X’Average Consolidation =’,F12.8///) 
FORMAT (4xX’KO =’ ,£15.6/4X’Beta a tal Sh al rd 
14X’Myv ee a JO Xie eke Sie Se 

14X’ GAMMA =’,F15.5/4xX’GAMMAW = Fass. 57% 

14X’U0 =! Fos, 5S: /4x? Vambaa "=" Fs .5/ 
14xX’Chi Se EV Tax NOT Or roots used =" 1S/ 
14X’No of elements ees / y 
FORMAT(13X’Excess pore pressure vs depth:’// 


113X’Depth’,&X’Depth’,4xX’NODE’,SX’Time t’/24X’Ratio’,10XG20.9 
17 /CSRF 11. SS, F107 4, 3X13, 424X615. 8) ) 

FORMAT(////5X’Error message from the *IMSLDPLIB for the final’ 
1,> “Tteratiron: “7 /sX’ TES IEJ2 IEJ3 PEY? IEY2‘°/3XS16) 


BREF RRKS SEEK SEK EK SETTER KEKE KTH KEK KKK KAAS KAKAKAHK ARETE EEE 
# IMSL Library # 

The following is a@ brief description of the Fortran 

routines used in this program. This information is 

available in the IMSL Library reference manual, edition 

6, 1980. The routines with single precision are stored 

in the public file *IMSLLIB, those with double precision 

are stored in *IMSLODPLIB. 

Routine name: MMBS JO 

Usage: FUNCTION MMBSJO(ARG,IER) 

Description of arguments: 

ARG - Input argument 

IER - Error parameter, terminal error 

Routine name: MMBSYN 


Usage: CALL MMBSYN(ARG,ORDER,N,YN, IER) 


Description of arguments: 


ARG - Input argument 

ORDER - Desired order, which must be greater than or 
equal to zere and less than one 

N - Number of results required, with which order 
Start from. “ORDER” to “ORDER S14 to ORDER S27) etic 

YN - Output vector of above results 

IER - Error parameter where 


‘129’: ORDER is out of range or ARG is less than 
the minimum allowable value 
“13206: HN VS Yess than or equal to zero 


“131’: YN overflow 
“132’: ARG is greater than the maximum allowable 
value 


* 
* 
* 
+ 
= 
* 
* 
* 
= 
* 
* 
* 
* 
+ 
3 
2 
* 
= 
= 
c 
+ 
x 
= 
* 
= 
* 
* 
* 
* 
* 
= 
* 
= 
=: 
* 
= 
= 
* 
* 
* 
= 
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KERR TK SESE ESET TE SH SSE SEA SAKE FHKE HK HHS EEK EREEATHK KEKE EE HE EAEE 


END 
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B.4 Edelmann's closed-form solution 


BESTS SST SSEC TSS SESS TEESE SETTERS ERE ESTAS eeeressesssaeesese 


ses EDELMANNS3 #28 


elasticity proportional to the Gepth 


analytical solution in 1853. The g0i11 
to drain at the top oniy. 


Input vaeriabdies: 
T © Time required, sec 
SIGMAO = & sma) initial effective str 


by the egn Myv=1/(C®Gammat2+S igmao 


GAMMA = Unit weight of soil, g/cmset3 
GAMMAW = Unit weight of water, g/cmts3 
Lal $o11 tnickness, cm 
K 
UO = Initial uniform pore pressure, g/ 
NOEL © WNumber of equal segments of the 
specified for the outputs 
VON(m) = Roots of eEGn: JO(V) YI(kKW)-YO 
where k 1S the variable LAMBD 
NORT = Total mumber of the above roots 


Use of tne program: 

a] First to compile the program by: 
RUN *®FORTGTEST SCARDS=EDELMANNS3 SPRIN 
2 Repeat the following for subsequent 
Remarks: 


Tne routines MMBSJO, MMBSJ1 and MMBSYN 


Back of th1i1s program 


ee ee ee) 


PESEKRSSL SETTERS TEST ETREL ESET ETE ETeA tt S 


IMPLICIT REAL®8(A-H,0-2) 


REAL*& LAMBDA,K,LIM,DSORT,OFLOAT,DEXP,DLO 
ENTEGERS DU GTA MON AIG eee eS a So Fg, 
INTEGER VF107)/°(S", 4:77, on 2 oe Sut 
INTEGER VF2(8)/'(SX:','T7HTN’,’eta °,’2:3K 
DIMENSION VN‘10!,21100) ,BUOVN‘(10),YOVN(10 
DIMENSION ALVYN( 10) ,2E(10),THETA‘10),2H(10 


BSPFFESSRESTrsrsesrsrsreerrsseserer srr srrszes estat 


= 
2 Reac alt input data 
2 


ERERCREC SRE SHR TRERSESERCL ER HYSET ST eESeERSEB EE 


READ(5,50)T,SI1GMAO,C,H,K 
READ(5,51)NOEL,NORT,UC,GAMMA, GAMMAW 
READ(S,S2)(VN‘1),]£1,NORT) 


Bezaeaeeeezreetreaseseseeezaereseeseezrseeesseeseeeezres 
= 

bd Variabie format input 

2 


VFI(S+:eDICITtnorRT) 
VF2(6)=D1G1T(NORT) 


Brrreeeerrersrsrerseersertererasrezerrrsesree xcs 


= 
2 
> Compute the average consolidation at t 
2 anag the % error at time Oo Tne thete 
s calculated 

z 

2 


EUuTSTRsevrrerereseresrsrressraerressesrrrseress 


ALPHS=S1GMA0/ (TC GAMMA) 
HLI=ZH/LLPHA+i DO 

LAMBDA=1 DO/DSORT( HAI! 
TT=C2K 2=GAMMAz7T/GAMMOW 


DHT=0.00 
CONV=O0 DO 


00 20 J=1,NORT7 
ALVN(1)=LAMBDA* VN(J) 

BJOVN: J] )=MMBSJO‘(VN()]),ITEJ1) 
BJUOLVN=MMBSJO‘(GLVN'1),1EJ2) 


BUIVN=EMMBSJI(VN( J), 1ES3) 
CALL MMESYNCVN( I.) ,60.00,2,YVN, TEY1?) 
CALL MMBESYN/ALVN(]),©0.00,1,YOLVWN,IEY2) 


YOVN(T)=YWN‘ TE 

YIVN=YVWVN(2) 
ZJOLVNZEBYOLVN-BUOVN( J} FYOLWN/YOWN']? 
ZI1VN=BUIVN-BYOVN())27YIVN/YOVN( 5 


* 


layer 


# One dimensional consolidation with the modulus of 


This program is to compute the excess pore pressure 
at various ceptnh at time t according to Edelmann’s 


1S assumed 


ess, g/cmts2 


) 


Tne coefficient of permeability, cm/sec 


cmts2 
soi) 


C = A dimensrontess constant, relates to compressibility 


where Z is the distance in cm from the surface 


profile 


(V) SI(KV)FO 
a2 wm the program 


useo 


T=-1 PAR=ID,SOURCE 
executions: 


DEBUG -LOAD#+=IMSLDOPLIB S=cGatatfile Groutput 


The % consolidation computed is in terms of 
im tne 
{ibrary wi11 be used Tneir cescriptions ere shown at the 


errers 


zuevage 


errere 


zaeenE 


sett 


Yement 


IMSL 
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8S Oy 
NOUC Kee eh 2) 
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THETA(JTIFO. 258DO2 ALVN(J)F®ALVYN( JT )27TT/ALPHA 
EXPT( II SDEXP'-THETACI)} 
BECTISAN( 11,2 EXPT) 
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DHT=DHT+AE(1)*ZOLVN 
CONV=CONV+AN(1)*ZOLVN/(2.D0*U0) 
20 CONTINUE 


LIM=DLOG(1.DO/LAMBDA)/2.D0 
ERR=(CONV-LIM)*100/LIM 


TEM=2.D0/(GAMMA*C) 
DELHF=TEM*UO#LIM#2.D0 
DELHT=DELHF - TEM*DHT 
AVCON=DELHT/DELHF 


RRR ESKER ETA AEE EERE ERK AKT KEKE KEKE EEKREKEEK EERE EERE ETE 


* * 
* Calculation of the excess pore pressure at various depth * 
* = * 


HRS KEK KEKE REESE ERE KEE KEE KEKE KKEKEREEEE EHR E REE EE EEEE 


NONOD=NOEL+1 
DIV=H/DFLOAT(NOEL ) 
Z(1)=0.D0 


DO 25 J=1,NONOD 
Z(S+1)=2(5)4DIV 
ZH(J)=Z2(5)/H 
ZA1=Z(J)/ALPHA+1.DO 
UE(J)=0.D0 


DO 25 L=1,NORT 
V=ALVN(L)*DSORT(ZA1) 
BUIVEMMBSJ1(V,1IEJ4) 
CALL MMBSYN(V,0.00,2,YV,IEY3) 
YIVEHYV(2) 
ZIVEBUIV-BJOVN(L)*YIV/YOVN(L) 
VE(J)=UE(J)+AE(L)*VFZ1V 

25 CONTINUE 


RHR KK SEHK ET KEK REHEAT HR EK KET AEH KEK AHKKS KE TK AK AA AKEKRKREKKAKKKKEKEEEE 


+ * 
* Output quantities * 
> = 


HRS HHH FHKE EKA EERE FSH EKEK KEKE SHAKE EHA KARE AKEAAARTE KEKE KE 


WRITE(6,60) 

WRITE(6,61) AVCON,ERR 

WRITE(6,65) K,C,SIGMAO,H,GAMMA, GAMMAW,UO, ALPHA, LAMBDA,NORT,NOEL 
WRITECG 62) T C2 (Cb), 2H Ll), b, UCL), l=1) NONCB) 

WRETECGY Gs hed, Lede nleus teva, LEV» Deve, LEYS 

WRITE(6,64) 

WRITE(6,VF1)(N,N=1,NORT) 

WRITE(6,VF2)(THETA(J),J=1,NORT) 


STOP 


SRF S SEEKS KAFFEE AEF KEE F SH KTS AEH EF KHK THK EA EA TEAK EK EAE KK KAKARKK KITE 


= = 
* Format statments z 
* * 


RRS EK EEK SEK EKER KSEE KA KEKE KEE K KEK SKA KES KHER AK AAA SEEHF HK RAT AEKFE EK KEE 


50 FORMAT(5G20.0) 
51 FORMAT(215,3620.0) 
S52 FORMAT(3F20.0) 


60 FORMAT(‘1’////4X’The One Dimensional Consolidation with “’, 
1’the Spatial Variation’/23X’of Compressibility’ 
1//13X’Analytical Solution by Edelmann (1853)'///) 

62 FORMAT(13X’Excess pore pressure vs depth: ’// 
113X’Depth’,6X‘Deptnh’ ,4X’NODE’,SX’Time t’/24X’Ratio’,10XG20.9 
V/ShSEKE UGS FIO 74 SKS, AXE 15. .8*)) 


61 FORMAT(4&X’Average Consolidation =’,F12.8 
tLe OO URCrOr Of ave JCOnSOll ate time O00 =, FiS cian Gauss fa) 

63 FORMAT(////SX’Error message from the *IMSLDPLIB for the final’ 
ioe tteratven. _/AS XX Le St IEJ2 be re ec} IESG LEY 1h 472 BEN siae: 


’ 


12X7(2X14)) 
64 FORMAT(///5X’The theta terms: ’) 
65 FORMAT(4&X‘K =" 24S .6/4xX%7C =e Fis.) 5¥ 


14X’°SIGMAO =’ ,F15.5/4&X’H COE tS. 57, 

14K ’GAMMA =’,F15.5/4X’GAMMAW SEs. S:/ 

14xX’U0 SC PAS eS /aexceip he 22 PF tSisS/ 
14X’Lambda =’,F15.5/4X’No of roots used =’,I5/ 
14X’No of elements em Iss //} 


HERSEK EEE KKH KSSH ASEH SKK AERA EHF EA SKA KEKE EK EHEKAH EK ATEEK HEHE 
# IMSL Library # 


The following is a brief description of the Fortran 
routines used in this program. This information is 
available in the IMSL Library reference manual, edition 
8, 1980. The routines with single precision are stored 
in the public file *IMSLLIB, those with double precision 
are stored in *IMSLDPLIB. 


Routine name: MMBS JO 
Usage: FUNCTION MMBSJO(ARG,IER) 
Description of arguments: 


ARG - Input argument 
IER - Error parameter, terminal error 


"+ # He He MH HH HHH HH HH HW HH HM OH 
HH HH HH HH HH HM HH HH HH HOH OM OH 


Routine name: MMBS J1 
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END 


Usage: FUNCTION MMBSJ1(ARG,IER) 
Description of arguments: 

ARG - Input argument 

IER - Error parameter, terminal error 
Routine name: MMBSYN 


Usage: CALL MMBSYN(ARG,ORDER,N,YN, IER) 


Description of arguments: 


ARG - Input argument 

ORDER - Desired order, which must be greater than or 
equal to zero and less than one 

N > Number of results required, with which order 
start from ‘ORDER’ to ‘ORDER+1’ to ‘ORDER+2’, 

YN - Output vector of above results 

IER - Error parameter where 


“129°: ORDER is out of range or ARG is iess than 


the minimum allowable value 
‘130’: N is less than or equal to zero 
°131’: YN overfiow 


‘132’: ARG is greater than the maximum allowable 


value 


RRSP SERA EKER TEKS AKER KAEK FREE A RES K THE HKEKEKEKE KEKE TEE KEKXKEXH KE 


etc 
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B.5 Datafile A.20: Eq. A4 
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BR ae a Me a eRe ok a Sk a aS a ea ak ok aE a Se iC ice ak ak ake eae ac ake oe ik Ok 8 ok ok OR a 


Varia 


To us 


+ © H# 4H © MH H H H He HH HH RH HH 


* 


alain lina ie eee ale Meee ee ee 2 a ae ee ee ee 


bies: CC 


DIVU 


E 
EO 
K 


Pp 
PO 
U 
XL 
e: RUN 


ie Fee et 


CMe feat] 


CE/AInt® 
load 


void ratio 


initial 


void ratio 


increment / number of material 


coefficient of permeability, k; 


or variable of 


"gamma*mv", 


where mv 


#This program is to prepare part of the data for A.20(Eq.A4) 


according to ADINAT users manual input section X.3 & X.4 


input points 


is the 


coefficient of volume compressibility 
effective stress 


initial 


effective stress 


excess pore pressure 


reduced 


depth, 


DEBUG -LOAD¥ 5=data 6=output 


IMPLICIT REAL*8(A-H,0-Z) 
REAL*8 DLOG,DLOG10,K(500) 
DIMENSION U(S00O), XL(500) 


A= 15D 
XL(1)=.1419D0 
XI) = "459200 
XL(3)=.8025D0 
XL(4)=1.1640D0 
XL(5)=1.5399D0 
XL(6)=1.9281D0 
XL(7)=2.3270D0 
XL(8)=2.7356D0 
XL(9)=3.1530D0 
XL(10)=3.5785D0 
XL(11)=4.0116D0 
XL(12)=4.4518D0 
XL(13)=4.8987D0 
XL(14)=5.3521D0 
XL(15)=5.8115D0 
XL(16)=6.2769D0 
KEC17) =6 "7478D0 
KIC 18) = 7225000 
XL(19)=7.70593D0 
XL(20)=8.1926D0 
CC=.4D0/DLOG(10.D0) 
DIVU=10.D0/20.D0 
U(1)=0.DO 
DOWOMT = 420 

UErI +s )=U0TT)+DIVU 
DO20T=1,20 
PO=.1D0+1.9D0*XL(I1) 
EO= 

EO1=1,.DO+EO 
DO30U=1,21 
P=P0+10.D0-U(J) 
KC UvIaCC ACR” EO 1) 
CONTINUE 
WRITE(6,60)I1.A 


WRITE(6,61)(U(LL).LL=1., 


SPO- F4D0* DEOGIO (107 50* PO) 


2) in CK CULE) erie= ati le) 


Wokhecer 61) CWO) etl =192 1) SOK) UL Sto) 


CONTINUE 
STOP 


FORMA TICS J T2E 10 20) 
EORMAUEG2 (BF 10. 1/7). SF 10). 11, SF 10.97/73 (8F 10.97 )) 


END 


computed from Eq. 
*FORTGTEST SCARDS=file SPRINT=-1 


6.3 


PAR=id,source,map 


in text 


2 ¢ % 4 UW 4 BU HM HB BH HM B- HB 


* 
* 
* 
* 


a i y wk sf an 


-)¢°6) S24: et 9.4 Oe 8 a 


crennbnnee 
, is Lik 7 7 
tna pa) jos. 10 
ee eg not 308 m4 
nv io ney 7 et, 
‘ented 1 ton Natrersn hedhaacn\ | 


1 ae ae ea Oa eerie 
ony at vin orerta "ve emma 


Pear mt €.a od Prey bazuqron , Fe 
an eptivor. oF Rae thal bak teas otis 3 


pt he 


nad i Ah ev! F rips 


SnVe Sea ee sola Sag sri ysoe 7) 
earths: evi taeyhe Taiz tril: 8 
| amen age Nady OH DH + 
BST cts an 
ry 2TaT l i 
SOND Ie uss ne 


tuaruore 


Llvasicassoatinacenuess sinh (es anna nw oe * 
Tam ae We nay 5 ae fetes 
)) fooe 3, or 90, 2040 8" 3A 
ie ee OEE (O08 du or 
ae a , coer. ft). 
* m3 ; ~ OOSOR> = 5. 
oe _ 9GatOs. E 
: oapaat if 
a 1. : odeeke e.! G 
Cae  gdvsae. tate 
ul probs rhe i 
i ee .  aeter, ‘eete), 
Rb  paceer Bee) 
2) vor i 
ee OREO bn te 
21 ia 
Ui] oorgea keer). can a 
f —oarsee. se( bt pdx: 
’ OCS B. Bel er PIR 
— Odeats, aeCar 
Ostet, ge Th), 
— OngaSS, THER ie . 
| resehe hee oe Ce 
OAeser. Be van 
| (oa. OF OO ID\OTS 33. . 
Oe a0. ane Oreuvid. 
Rye ah Bp 
uiaae . OS tat tor . 
he pass ty 
‘wie, (i, Se re] 
+ it) deedoe-t 
JApaeo. or 101 OSES 


04oe 


ih nS |: 
' at Pa a 
ns 


. AS ‘ae te. ere ti 
ATR, erin te ht 


OOO) O11: OO O1OLGiG)' OC OG) GG) 


10 


30 


20 


60 
61 
62 
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BERG RR IRR ee SEARO ROR SRS RCE SC PC aKa SR SEK 6 CRE Re SEE OR RT oe TN RSE Re See EL we 
* * 
* #THis program is to prepare part of the data for AS ZO GE CAB) z 
* according to ADINAT users manual input section xX.3 & X.4 x 
* Variables: CC =" Cc/lini© * 
- DIVU = load increment / number of material input points * 
2 E = void ratio Ps 
~ EO = “itettade Vow daratso 2 
* K = variable of "k*(1+eo) / i+e" % 
* MV = variable of "gamma*mv*(i+e) / 1+e0" yy 
* P = effective stress 3 
* PO = initial effective stress = 
* U = excess pore pressure 
* XL = reduced depth, computed from Eq. 6.3 in text * 
* To use: RUN *FORTGTEST SCARDS=file SPRINT=-1 PAR=id,source,map ar 
* * 

* 


Mee eo kee eae ao ak ke ok kk hk kk OOK Ke ee OR KE 


DEBUG -LOAD¥ S5=data 6=output 
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B.6 Datafile A.20: Eq. A8 


IMPLICIT REAL*8(A-H,0-Z) 
REAL*8 DLOG,DLOG10,.K(500) ,MV(500) 
DIMENSION U(500) ,XxL(500) 


A=1.D0 
XL(1)=.1419D0 
XL(2)=.4592D0 
XL(3)=.8025D0 
XL(4)=1.1640D0 
XL(5)=1.53993D0 
XL(6)=1.9281D0 
XL(7)=2.3270D0 
XL(8)=2.7356D0 
XL(9)=3.1530D0 
XL(10)=3.5785D0 
KE Gini 40 6DO 
XL(12)=4.4518D0 
XL(13)=4.8987D0 
XL(14)=5.3521D0 
KEES) =5. 84115D0 
XL(16)=6.27693D0 
XL(17)=6.7478D0 
¥iE618 )=7..2250D0 
XL(198)=7.7058D0 
XL(20)=8.1926D0 
CC=.4D0/DLOG(10.D0) 
DIVU=10.D0/20.D0 
U(1)=0.D0 

DOA ONIL=1 <20 


UGII+1)=U(TE)+DIVU 


POZOL=dy, 


20 


PO=.1D0+1.9D0*xL (1) 
EO=.9D0- .4D0*DLOG10( 10*PO) 
EO1=1.D0+EO 


DO30uU=1 
P=PO+10. 


oa 


DOLCE 


E=.9D0- .4D0*DL0G10(10*P) 


E1=1.DO+ 
MV(J)=CC 


Ki pECCSEOA/(P*E1*E4)) 


CONTINUE 
WRITE(6, 


Writer G. 64) (UCLE) LU 4. 


E 
/(P*EO1) 


60)1,A 


2 AICS) USS.) 21) 


Write C62 Ge \(UICEE) sy bE= I 23) CMMCUe i Ean.) 


CONTINUE 
STOP 


FORMAT(3X,12,F10.0) 


FORMAT(2(8F10.1/),5F10.1,3F10.9/3(8F10.9/)) 
FORMATC2CGrRHOe 14) bk 1On Teer 1OnS/ Ser iOno7)e) 
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B.7 Datafile A.20: Eq. A9 


SE ak kk kok oo ok ok ok ook OF Fo rk Kk OO OR OO OK 


DEBUG -LOAD4 5=data 6=output 


* 


WO oe ok ok OK ok kak ok kk ok Oo ok ok ok OF ok ok ok kk ko ok kook 


IMPLICIT REAL *8(A-H,0-Z) 
REAL*8 DLOG,DLOG10O,K(500) ,.MV(500) 
DIMENSION U(500) 


A=1.DO 
DIVU=10.D0/20.D0 
U(1)=0.DO 


DO1OI1=1,20 

UC II+41)=U(II)+DIVU 

CC=.4D0/DLOG(10.D0) 

DIVZ=8 .43786D0/20.D0 

Z=DIV Zeb oO 

DO201=1,20 

DO30U=1,21 

P=10.1D0+1.9D0*Z-U(UJ) 

E=.9D0- .4D0*DL0G10(10.D0*P) 

E1=1.DO+E 

MV(JU)=CC/P 

K(J)=MV(J)/(E1*E1) 

CONTINUE 

WRITE(6,60)I1.A 

WRpme Ger CepnCUiC sia 2 4) CK CU) uw 24) 
WRimMeGGe Gi rGUiGnL) El =1. 94). (MY CUC). JUu=te 2s) 
Z=Z+DIVZ 

CONTINUE 

STOP 

FORMAT(3X,12,F10.0) 
FORMAT(2(8F10.1/),.5F10.1,3F10.9/3(8F10.8/)) 
FORMAT(2(8F10.1/),5F10.1,3F10.9/3(8F10.9/)) 
END 


To use: RUN *FORTGTEST SCARDS=file SPRINT=-1 


interpolation points 


PAR=id,source,map 


* 

* #THiS program is to prepare part of the data for A.20(Eq.A9) 
a according to ADINAT users manual input section X.3 & X.4 
* Variables: FCC = Cofinie 

ba DIVU = load increment / number of 

a DIVZ = reduced layer depth / number of elements 
* E = void ratio, e 

* K = Vvartablemofiwlarse +e" 

* MV = variable of "gamma*mv*(i+e)" 

* Pp = effective stress 

* U = excess pore pressure 

* 74 = reduced coordinate, ‘’b’ 

x 

* 
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B.8 Datafile A.21: Eq. A4 


Oe Ok Ok OO OO OK 
* * 
* #THiS program is to prepare part of the data for A.21(Eq.A4) z 
* according to ADINAT users manual input section X.3 & X.4 bs 
+ Variables: DIM = load increment / number of material input points * 
> E = void ratio ¥ 
x EO = initial void ratio % 
* K = k ( = gamma*mv*Cv ) % 
* MV = variable of "gamma*mv" t 
ny P = effective stress % 
* POO = variable of "-lambda * initial effective stress" * 
* U = excess pore pressure 3 
nah = reduced depth: Eq.31 from Gibson etal (1981) = 
* To use: RUN *FORTGTEST SCARDS=file SPRINT=-1 PAR=id,source,map * 
* DEBUG -LOAD4 5=data 6=output * 
* * 
oe ok ok oe ok ook ke oko ok ok kook ok aka oka oko ko kok ok ok ok ok ok Fo ok ok ok kok ok ok ok Ok Ok OF Ok OOK 


IMPLICIT REAL*8(A-H,0-Z) 

REAL*8 DEXP,K(500) ,MV(500) ,XL(500) 
DIMENSION U(500) 

A= OO 

XL(1)=.0521D0 

XL(2)=.1587D0 

XL(3)=.2684D0 

XL(4)=.3813D0 

XL(5)=.4974D0 

XL(6)=.6168D0 

XL(7)=.73986D0O 

XL(8)=.8658D0 

XL(9)=.995D0 

XL(10)=1.128D0 

KU Gt) a44265D0 

KL(12)=1.405D0 

XL(13)=1.548D0 

XL(14)=1.695D0 

XL(15)=1.845DC 

XL(16)=1.998D0 

XL (17) =2. 155D0 

XUG8) S243 15D0 

Xi CNS 47 DS 

XL(20)=2.642D0 

DIVU=10.D0/20.D0O 

U(1)=0.DO 

DO1OII=1,20 

U(I1+1)=U(I1)+DIVU 

DOZOL=1, 2O 

POO=-.332D0*1.74D0*XL(1) 
EO=1.49D0+2.34D0*DEXP(POO) 
EO1=1.DO+EO 

DO3OU=1,21 

P=- 2362002 (4 74 D072 xb Gl) 10, DO-UiCU)) 
E=1.49D0+2.34D0”*DEXP(P) 
AM=..332D0*(E-1.49D0) 

MV (JU) =AM/EO1 

K(J)=MV(JU)* .001944D0 

CONTINUE 

WRITE(E.60)1,A 

WRinkE CGE G1) (UC Isis Wa ot tae CK (we ie Ul =o) 1) 
WeieneGGr Goo GU (EEO, biadin® te OMV GUE Dm Ul =e) 
CONTINUE 

STOP 

FORMAT(3X,12,F10.0) 

EORMAMCS (SEdO<H/ ) 45F AO od. SEAM 5S /S0SEta 467 ))) 
FORMAT(2(8F10.1/),5F10.1,3F10.9/3(8F10.9/)) 
END 


“peente dvitsedte Tarntat * shamets* “to bea — OF 


(Rr) late onac!d sort spk a Berwees + IK ? 
oon, go UE DreHAG T- <n ‘ wy rt betnnni tos : 
a . Par wk P 4) i 


. POPTTPeTTL Tt er . 
4 r 


ix te q an 


Me ctemes _ a pet cM os 
le 


ih t 


a 6 - « 
5) —. a 
a 1 
a 


4 {6 : 


v 


‘ , oh 
ae = os 
: . = 
eeenve, (anes BR TTSTIaMS 
adel ve (¢ nae a a ~ Rad 


ovtaer tL pUntes | 
re 5 ; 


ti vux-ooet oat oad 
(005 an TaOdeE 


Pyey¥-od. ot +t 1) 227 OGRe. Hee 
19) TKID"OOKE & 


Sys Aik (rs. te 
Recor ald re; TEP) 


GOO OOO OOOO Oa OG Ooi 


30 


60 


62 


BL? Datarliic A.21 > gq, As 
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a ee eee eee eee eee eee ee ee ee ee a ee ee a ee a 


Variables: DIVU 


t+ 4 © © 2 £2 BOM HF BR DR HH Ht 


To use: RUN *FORTGTEST SCARDS=file SPRINT=-1] 
DEBUG -LOAD¥ 5=data 6=output 


* 


REE EEK a EEE KK 


IMPELC IR REARS GASHMO=Z ) 
REAL*8 DEXP,K(500) ,MV(500) ,XL(500) 
DIMENSION U(500) 


XL(1)=.0521D0 
XL(2)=.1587D0 
XL(3)=.2684D0 
XL(4)=.3813D0 


XL(5)=.4974D0 
XL(6)=.6168D0 
XL(7)=.7396D0 
XL(8)=.8658D0 


XL(9)=0.985D0 

XL(10)=1.128D0 

XE(14)=42265D0 

XL(12)=1.405D0 

XL(13)=1.548D0 

XL(14)=1.695D0 

XL(15)=1.84550 

XL(16)=1.9898D0 

XC 7 =2 455D0 

XL(18)=2.315D0 

XL(19)=2.477D0 

XL(20)=2.642D0 

DIVU=1C.D0/20.D0 

U(1)=0.DO 

DOIOII=1,20 

Vem 1 =U LI) +D1VvuU 

DO20I=1,20 

POO==.332D0* 1. 74007 XE (1) 
EO=1.49D0+2.34D0*DEXP( POO) 

EOd=1-DOTEO 

DO30U=1,21 

P= —) 33200 (4a 406k XE OL) 1, DO=UICG 
E=1.49D0+2.34D0*DEXP(P) 

£1=1.DO+E 

AM=.332D0*(E-1.498D0) 

MV(JU)=AM/EO1 

K(JU) =AM* .00194400*E01/(E1*E1) 

CONTINUE 

WRITE(6,60)1.,A 

eager Metts (Sap OKC E NUL Soh iy) EPG) wks} sah, 
WRT Gun Gu) (rl) C (ates) ie eile —ealare a) ea MVC) (ee) ee cli ieeeat) 
CONTINUE 

STOP 

FORMAT (SX) l2eF 10203 

FORMAT ( 2(8F 10. 1/)).5F 10.1, SE IM 6/3085 14-767)) 
FORMAT (COCBE 10. 1/7.) SPIO st SSF 1049/3 08F 10297 ))) 
END 


#THis program is to prepare part of the data for A.21(Eq.A8) 
according to ADINAT users manual input section X.3 & X.4 
load increment / number of material input points 


E = void ratio 

EO Si nitratevoldyra tio 

K = variable of "k*(1+eo) / 1+e" 

MV = variable of "gamma*mv*(i+e) /(1+e0)” 

Pp = effective stress 

POO = variable of "-lambda * initial effective stress" 
U = excess pore pressure 

XL = reduced depth: Eq.31 from Gibson etal (1981) 


PAR=id,source,map 


* 
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B.10 Datafile A.21: Eq. A9 


We eee oo ok ok ok ok ok a ok ok ok ok ok ok kok ok ok ok oo ak ok ok kk ok kk kk kk ek ok 
* 


DEBUG -LOAD¥ 5=data 6=output 


* 


* 

* 
2 # This program is to prepare part of the data for A.21(Eq.AQ9) ‘3 
: according to ADINAT users manual input section X.3 & X.4 of 
* Variables: DIVU = load increment / number of interpolation points * 
be DIVZ = reduced layer depth / number of elements = 
= E = void ratio, e a 
* K = variable of "k / 1+e" . 
* MV = variable of "gamma*mv*(it+e)" od 
* P = effective stress 7 
* U = excess pore pressure bs 
* Z =e Reaouced sCOonGinatenus on a 
* To use: RUN *FORTGTEST SCARDS=file SPRINT=-1 PAR=id,source,map * 
* * 


We ook kk ek ok ok ok ek ok kk ok kok ok ok ko ok OR ak ok ok OF ok oe ok ok ok kk ok oko 


TM PSEC RE Aas 8) GAS OZ.) 
REAL*8 DEXP,K(500) ,MV(500) 
DIMENSION U(500) 


A=1.D0 
DIVU=10.D0/20.D0 
U(1}=0.D0 


DOWOiT=40 2O 
U(II+14)=U(II)+DIVU 

DIVZ=2_728D0/20.DO 

Z=DIV2/2.D0 

DO2ZOI=1.20 

DO30U=1,21 

P=-.332D0*(1.74D0*2Z+10.DO0-U(U)) 
E=1.49D0+2.34D0*DEXP(P) 

E1=1.DO+E 

MV(JU)=.332D0*(E-1.49D0) 
K(JU)=MV(J) * .001944D0/(E1*E1) 

CONTINUE 

WRITE(6,60)1,A 

Meno SW) CUe) Ss) Sik, Ray, hts ah, apy) 
WRIME(C662) (UCR) hG= 124). (MVC) SUE ts 24) 
Z=Z+DIVZ 

CONTINUE 

STOP 

FORMAT(3X,12,F10.0) 
FORMAT(2(8F10.2/).5F10.2,3E11.6/3(8E11.6/)) 
FORMAT(2(8F10.2/),5F10.2,3F10.9/3(8F10.9/)) 
END 
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